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ABSTRACT 

We present a detailed analysis of a large sample of 31 low-redshift, mostly radio-quiet type 1 quasi-stellar objects (QSOs) observed 
with integral field spectroscopy to study their extended emission-line regions (EELRs). We focus on the ionisation state of the gas, 
size and luminosity of extended narrow line regions (ENLRs), which corresponds to those parts of the EELR dominated by ionisation 
from the QSO, as well as the kinematics of the ionised gas. We detect EELRs around 19 of our 31 QSOs (61%) after deblending 
the unresolved QSO emission and the extended host galaxy light in the integral field data with a new dedicated algorithm. Based 
on standard emission-line diagnostics we identify 13 EELRs to be entirely ionised by the QSO radiation, 3 EELRs are composed of 
H ii regions and 3 EELRs display signatures of both ionisation mechanisms at different locations. The typical size of the ENLR is 
~10kpc at a median nuclear [Om] luminosity of log(L([0 m])/[ergs -1 ]) = 42.7 ± 0.15. We show that the ENLR sizes are least a 
factor of ~2 larger than determined with the Hubble Space Telescope, but are consistent with those of recently reported type 2 QSOs 
at matching [Om] luminosities. The ENLR of type 1 and type 2 QSOs therefore appear to follow the same size-luminosity relation. 
Furthermore, we show for the first time that the ENLR size is much better correlated with the QSO continuum luminosity than with 
the total/nuclear [O m] luminosity. We show that ENLR luminosity and radio luminosity are correlated, and argue that radio jets even 
in radio-quiet QSOs are important for shaping the properties of the ENLR. Strikingly, the kinematics of the ionised gas is quiescent 
and likely gravitationally driven in the majority of cases and we find only 3 objects with radial gas velocities exceeding > 400kms~' 
in specific regions of the EELR that can be associate with radio jets. In general, these are significantly lower outflow velocities and 
detection rates compared to starburst galaxies or radio-loud QSOs. This represent a challenge for some theoretical feedback models 
in which luminous QSOs are expected to radiatively drive an outflow out to scales of the entire host galaxy. 

Key words. Galaxies: active - quasars: emission-lines - Galaxies: ISM 



1. Introduction 



Extended warm ionised gas with sizes of several tens of kpc 
— i around luminous radio-loud quasi-stellar objects (QSOs) was 
O^l initially discovered by W ampler et al.| ( |1975| l and [Stockton 
(|1976|l. Fo llow-up longslit spectroscopy presented by Boro- 
|son & Oke] ([1984} and |Boroson et al.| ( |1985| ) confirmed that a 
• *h substantial part of the extended ionised gas is predominantly 
f\ photoionised by the hard radiation of the active galactic nu- 



cleus (AGN) as indicated by the similar line ratios compared 
to the compact narrow-line region (NLR). |Stockton & MacK^] 
|enty| ( |1987) conducted a large ground-based narrow-band imag- 
ing survey of 47 luminous QSOs at z < 0.5 in the [O in] /15007 
line (shortly [Om] hereafter) and detected ionised gas beyond 
10 kpc around 25% of the objects, predominantly radio-loud 
QSOs. These nebulae are often referred to as extended emis- 



* Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck-Institut fiir Astronomie and the Instituto de Astroffsica de An- 
dalucfa (CSIC) 

** Table 3, 4 and 7 are available in electronic form at 
http ://w w w. aanda. org 



sion line regions (EELRs) when their [O ra] line luminosity ex- 
ceeds 5 x 10 41 ergs -1 . A strong alignment effect between the 
major EELR axis and the radio axis was found for RLQs at 
Z > 0.6 ( |McCarthy et al.|1987 1, suggesting that radio jets have a 
strong impact on the EELR properties and may directly enhance 
their brightnesses. Furthermore, EELRs are preferentially found 
around radio-loud QSOs with a systematically lower broad line 
region (BLR) metallicity (Fu & Stockton 2007). Thus, the nature 
of EELRs and the origin of the gas are not yet clear and remain 
strongly debated. 

The properties of the extended ionised gas around radio- 
quiet QSOs has been studied much less intensively. Bennert 
|et al.| ( p"002) > imaged a sample of 7 objects in the [Om] line 
with the Hubble Space Telescope (HST) and detected ionised 
gas on scales of up to 10 kpc, significantly smaller than the 
EELR around radio-loud QSOs. Including the NLR sizes of low- 
luminosity Seyfert galaxies, |Bennert et al.| ( |2002| ) established 
a correlation between the [O in] luminosity and the size of the 
NLR, often referred to as the extended NLR (ENLR) when its 
size exceeds 1 kpc. The adopted distinction between the terms 
ENLR and EELR throughout the paper is that an (E)NLR des- 
ignates only the AGN photoionised phase of the emission-line 
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region, whereas an EELR refer to the entire emission-line re- 
gion, independently of the ionisation mechanism. Another large 
HST narrow-band snapshot survey targeting Seyfert galaxies 
confirmed the existence of such a relation (Schmitt et al. 2003b ), 
but the exact slope and zero-point still remain controversial. Po- 
tential caveats of the HST images are the low sensitivity and 
the lack of spectroscopic information to verify the ionisation 
mechanism of the gas to be AGN photoionization. For example, 



longslit spectroscopy of Seyfert galaxies allowed Bennert et al 



(2006a b) to estimate a characteristic distance of 0.7-4 kpc from 
the nucleus at which the photoionization of star forming region 
started to dominate against the AGN photoionisation. 

The study of extended ionised gas around luminous QSOs 
is significantly hampered by the severe contamination from the 
emission of the bright nucleus. Recent studies therefore focused 
on rare type 2 QSOs for which the emission of the nucleus is 
thought to be obscured by an optically thick torus following 
in the Unified Model of AGN ( |Antonucci||1993| l. For exam- 
ple, |W^MStfalt^(|20T0| detected a huge EELR (180 kpc) 
around a luminous type 2 QSO (z = 0.399) and determined an 
ENLR size of ~25 kpc from longslit spectra using emission line 
diagnostics. Furthermore, small samples of type 2 QSOs were 
recently studied by means of longslit (Gree neet al.|201 ljfV illar- 



Martin et al.|2011b] l and integral field spectroscopy ( Humphrey 



et al. 20 1 0]> . Typical ENLR sizes of ~ 10 kpc were reported, 



almost 2 times larger than the ENLRs of unobscured (type 1) 
QSOs observed with HST. 

The extended warm ionised gas around AGN has attracted 
further attention because it may allow to directly study the im- 
pact of an AGN on its surrounding host galaxy. Complex kine- 
matics in the EELR are found around radio-loud QSOs and ra- 
dio galaxies with line widths up to 1000 km s and high velocity 



clouds that exceed the escape velocity of the galaxy (e.g. Fu & 
Stockton||2009l |Holt et al.||2006| |2003) [McCarthy et al.||1996 



Nesvadba et al. 2006; Tadhu nter et al.|2001||Wl ar-Martm etaL] 
1999). These QSOs are therefore expected to expel a large frac- 



tion of gas from the host galaxy. This is usually interpreted as a 
clear signature for AGN feedback which has been invoked to ex- 
plain the observed scaling relation between the black hole (BH) 
mass and host galaxy properties (e.g. Gebhardt et al. 2000, Giil- 
tekin et al.||2009[ |Haring & Rix||2004[ |Magorrian et al.||1998[ 



Tremaine et al.| 2002 1, and to efficiently quench star formation. 



This is required to prevent the excess of massive young galaxies 
that is often seen in cosmological simulations (e.g. |Croton et al. 
2006 ; Khalatyan et al. 2008 1. Given the short duty cycle of radio- 
loud AGN and the low radio-loud fraction of the order of 10%, 
it remains still unclear whether the radio-loud mode of AGN ac- 
tivity can be the primary channel for AGN feedback. Further- 
more, the coupling of kinetic jet energy to the ambient medium 
may be significantly less efficient than currently required by nu- 
merical simulation ( Holt et aT7||201 1] >. Extremely complex ex- 
tended gas kinematics similar to those in radio-loud counter- 
parts have also been confirmed in individual radio-quiet QSOs 
like the ultra luminous infrared galaxy (ULIRG) Mrk 23 1 jAalto 
et al.|)2012l|Hamilton & Keel||1987[ |Lipari et al.|[2009l |Rupke 



& Veilleux|2011| l, but the majority of radio-quiet type 2 QSOs 



appears to have more quiescent gas kinematics ( [Greene et al.| 
|2011 1 [Humphrey et al.||2010| . An exception are double-peaked 
[O in] AGN, which were initially thought to be binary AGN can- 
didates, but Fu et al. (2012 1 reported that the majority of them 
are simply AGN with a luminous ENLR based on integral-field 
spectroscopy. The line splitting of several hundred to thousand 
km s _1 suggests the presence of either significant outflows driven 
by the AGN or severe inflows as a result of merging processes. 



Anyway, the fraction of double-peaked [O m] AGN is < 1 % ( |Liu| 
et al.|2010[ [Smith et al.|2010l|Wang et al.|2009 1. 

Given that type 1 QSOs are much more frequent and allow to 
infer basic AGN properties from the spectrum of the nucleus, we 
took the challenge to deal with the deblending of the QSO and 



extended host galaxy emission. In our early work (Husemann 
|et al.|200"8| hereafter Hu08), we studied a sample of 20 mainly 
radio-quiet QSOs observed with an integral field unit (IFU) and 
found that 8 of the 20 QSOs are surrounded by an EELR. We em- 
phasised that its presence is linked to the spectral properties of 
the nucleus itself, in particular the equivalent width (EW) of the 
Fe ii emission and the width of the broad HB line. Since then, we 
further extended our sample with additional observations of 1 1 
radio-quiet QSOs. In this paper we present a detailed analysis of 
the EELR ionisation mechanisms, the ENLR size and luminosity 
as well as the ionised gas kinematics for our entire sample. 

The paper is organised as follows. In section [2] we report 
on the characteristics of the sample including its selection and 
host galaxy properties based on broad-band imaging. The IFU 
observations, data reduction and QSO-host deblending are de- 
scribed in Section [3] Basic properties of the AGN are deter- 
mined from the QSO spectra in Section|4]and we apply standard 
emission-line diagnostics to infer the ionisation mechanism of 
the extended gas in Section [5] We describe the properties of the 
ENLR in Section [6] and discuss them in comparison to previous 
results in the literature. Finally, we present the kinematics of the 
extended gas in Section|7] 

Throughout the paper we assume a standard ACDM cos- 
mological model with Hq = 70kms~ 1 Mpc, Q m = 0.3, and 
Q A = 0.7. 



2. The QSO sample and characteristics 

2.1. Sample selection 

Our QSO sample consists of 31 objects in the redshift range 
0.06 < z < 0.3, which we observed with the Potsdam Multi- 
Aperture Spectrophotometer (PMAS, |Rofh et al.|2005) > mounted 
to the Cassegrain focus of the 3.5 m telescope at the Calar Alto 
observatory. The QSOs were drawn from two samples that 
were constructed with different selection criteria. An initial 
sample of 21 low-redshift QSOs was randomly selected from 



the Palomar-Green Survey (PG, [Schmidt & Green||1983|l and 
from the Hamburg/ESO survey (HES, |Wisotzki et al.||2000) 
with an absolute magnitude of Mb < -23 (computed with 
Hq = 50kms _1 Mpc -1 ), a commonly used yet arbitrary divi- 
sion line between Seyfert galaxies and QSOs, and a declination 
of 6 > -10° to be observable from the Calar Alto observatory. 
An additional set of 10 QSOs were drawn from the Sloan Dig- 
ital Sky Survey (SDSS) DR6 ( |Adelman-McCarthy et alT2 008). 
These QSOs were selected to be brighter than g < 17 within the 
redshift range 0.13 < z < 0.3 and to have a [Om]/H/? flux den- 
sity ratio at their peaks greater than 1.5 as measured from the 
SDSS spectra. This last criterion is based on early results from 
the initial sample (Hu08) and should give preference to QSOs 
surrounded by extended ionised gas. The complete list of our 
observed objects is given in Table [T] including information on 
some of their main characteristics. 

The redshift-luminosity diagram of our observed QSO sam- 
ple is shown in Fig.[T] We computed absolute magnitudes (Mb) 
from the apparent B band magnitudes taking into account Galac- 
tic foreground extin ction (|Schlegel et al.|1998) > and an appropri- 
ate ^-correction (Wisotzki 2000|l. The B band magnitudes and 
redshifts were either taken from the available catalogues or di- 
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Table 1. Characteristics of the QSO sample. 



Name 


Alt. Name 


a (J2000) 


S (J2000) 




m B b 




/l.4GHz d 


fs GHz 6 


R ! 


Class. 8 


PG 0026+129 




OO h 29 m 1358 


13°16'05'.'0 


0.142 


14.95 


-24.24 


7.1 


5.1 


1.08 


RI 


PG 0050+124 


I Zw 1 


OO h 53 m 3459 


12°41'35'.'0 


0.061 


14.39 


-22.86 


8.3 


2.6 


0.33 


RO 


PG 0052+251 




00 h 54 m 52?2 


25°25'39'.'0 


0.155 


15.42 


-23.97 


7.1 


0.7 


0.24 


RO 


SDSS J0057+1446 




00 h 57 m 09?9 


14°46'10'.'2 


0.172 


15.58 


-24.05 


2.4 


1.3+ 


0.48 


RO 


HE 0132-0441 


RBS 219 


01 h 35 m 27?0 


-04°26'36'.'2 


0.154 


16.25 


-23.12 


7.4 


2.9 


2.03 


Rl 


SDSS J0155-0857 




01 h 55 m 30;0 


-08°57'04'.'l 


0.165 


16.83 


-22.71 


< 0.9 


< 0.5* 


<0.58 


RO 


HE 0157-0406 




01 h 59 m 4950 


-03°52'00'.'4 


0.218 


17.09 


-23.11 


< 0.9 


< 0.5* 


<0.75 


RO 


PG 0157+001 


Mrk 1014 


01 h 59 m 50?2 


00°23'41'.'0 


0.164 


15.20 


-24.32 


26.2 


8.0 


2.12 


RI 


SDSS J0836+4426 




08 h 36 m 5859 


44°26'02'.'3 


0.254 


15.54 


-25.03 


10.0 


5.3 1 ' 


1.93 


RI 


SDSS J0948+4335 




09 h 48 m 59?5 


43°35'19'.'0 


0.226 


16.81 


-23.48 


< 0.9 


< 0.5 + 


<0.56 


RO 


SDSS Jl 131+2632 




ll h 31 m 09?2 


26°32'07'.'9 


0.244 


17.41 


-23.06 


2.7 


1.4* 


2.90 


RI 


HE 1217+0220 


PKS 1217+023 


12 h 20 m ll?9 


02°03'41'.'7 


0.239 


17.23 


-23.20 


652.0 


440.0 


757.18 


RL 


HE 1226+0219 


3C 273 


12 h 29 m 06?7 


02°03'07'.'9 


0.158 


11.42 


-28.01 


55296.0 


37000.0 


302.22 


RL 


SDSS J1230+6621 




12 h 30 m 2252 


66 21'54'.'7 


0.184 


16.75 


-23.05 


< 2.5 


< 1.3 1 " 


<1.46 


RO 7 


HE 1228+0131 




12 h 30 m 50 ! ;0 


01 o 15'21'.'8 


0.117 


14.82 


-23.92 


1.4 


0.2 


0.04 


RO 


SDSS J1230+1 100 




12 h 30 m 54?l 


H o 00'll'.'2 


0.236 


16.61 


-23.78 


< 1.0 


< 0.5 1 


<0.52 


RO 


PG 1427+480 




14 h 29 m 43;l 


47°47'27'.'0 


0.221 


16.33 


-23.90 


< 0.9 


0.0 


0.02 


RO 


SDSS J1444+0633 




14 h 44 m 14s7 


06°33'06'.'8 


0.208 


17.10 


-22.99 


< 1.0 


< 0.5 1 


<0.81 


RO 


HE 1453-0303 




14 h 56 m ll?5 


-03°15'27'.'8 


0.206 


15.68 


-24.39 


2.7 


1.4* 


0.58 


RO 


PG 1545+210 


PKS 1545+210 


15 h 47 m 4355 


20°52'17'.'0 


0.266 


16.05 


-24.63 


1560.0 


720.0 


418.15 


RL 


PG 1612+261 


TON 256 


16 h 14 m 13?2 


26°04'16'.'0 


0.131 


16.00 


-23.00 


17.9 


5.1 


2.81 


RI 


SDSS J1655+2146 




16 h 55 m 51?4 


21°46'01'.'8 


0.154 


16.00 


-23.37 


6.6 


3.5* 


1.94 


RI 


PG 1700+518 




17 h 01 m 2459 


51°49'21'.'0 


0.292 


15.43 


-25.48 


21.6 


7.2 


2.36 


RI 


PG 2130+099 


Mrk 1513 


21 h 32 m 27*8 


10°08'19'.'0 


0.061 


14.62 


-22.63 


6.0 


2.0 


0.32 


RQ 


HE 2152-0936 


PHL 1811 


21 h 55 m 0L5 


-09°22'24'.'5 


0.192 


14.26 


-25.63 


1.2 


0.6* 


0.07 


RQ 


HE 2158-0107 




22 h 01 m 03!l 


-00°53'01'.'0 


0.213 


16.69 


-23.45 


1.6 


0.8* 


0.88 


RQ 


HE2158+0115 


RBS 1812 


22 h 01 m 29?7 


01°29'56'.'l 


0.160 


16.63 


-22.83 


< 0.9 


< 0.5* 


<0.48 


RQ 


PG 2214+139 


Mrk 304 


22 h 17 m ll?6 


14°14'28'.'0 


0.067 


14.98 


-22.48 


< 2.5 


0.2 


0.05 


RQ 


HE 2307-0254 




23 h 10 m 29?7 


-02°38'13'.'0 


0.221 


16.90 


-23.34 


< 2.5 


< 1.3* 


<1.68 


RQ? 


HE 2349-0125 


PKS 2349-014 


23 h 51 m 56?0 


-01°09'13'.'6 


0.174 


16.03 


-23.63 


1623.2 


622.7 


355.46 


RL 


HE 2353-0420 




23 h 56 m 32?9 


-04°04'00'.'7 


0.229 


17.23 


-23.09 


< 2.5 


< 1.3* 


<2.27 


RQ? 



Notes. w QSO redshift taken from the PG, HES, or SDSS catalogues. (b) Apparent B band magnitudes taken from the catalogues or sy nthesised 
from the SDSS spetrum. (c) Absolute magnitudes M B taking into account Galactic foreground extinction and an appropriate fc-correction ( Wisotzki 
2000 [ l. <d) Radio fluxes at 1.4 GHz in mjy as measured by the NVSS all-sky survey or the FIRST survey when available. <c) Radio fluxes at 
5 GHz in mjy taken from Kellerman n et al.| (JT989 ) and an unpublished catalog (Gopal-Krishna et al. private communication), if available. 
Values marked with f were estimated from the radio flux at 1.4 GHz assuming a spectral index of a = -0.5. (f) Ratio R of the 5 GHz radio to 
optical flux density at 4400A. Th e optical flux density / V (4400A) in mjy was estimated from the opti cal B band magnitude w ith the relation 
m B = -2.5 log / y (4400 A) - 48.36 ( jSchmidt & Green|l983| . (g) Radio classification of the QSO following [Kellermann et aT]|l989| into radio loud 
(RL) and radio quiet (RQ) and radio intermediate (RI) QSOs. An unambigous classification for a few object cannot be made due to upper limits 
in R. 




Fig. 1. Absolute B band magnitude against the redshift for all objects 
in our QSO sample. 



rectly measured from publicly available SDSS spectra. The ab- 
solute magnitudes are in the range of -26 < Mb < -22.3 with 
a median value of -23.45. The median redshift of the sample 
is z = 0.184. Thus, the majority of our objects are luminous 
low-redshift QSOs. 



2.2. Radio properties 

Radio-loud and radio-quiet QSOs (RLQs and RQQs, respec- 
tively) are often distinguished based on the R parameter, which 
is defined as the ratio of the radio (at 6 cm/5 GHz) to the optical 



flux density (at 4400 A). We adopt the classification of Keller- 
|mann et aL] ( |1989| l, who separated RQQs with R < 1, RLQs 
with R > 10 and a few radio-intermediate objects in between 
the two classes. Not every object of our sample has been ob- 
served at 5 GHz, but the NRAO VLA Sky Survey (NVSS [Con 



don et al. 1998 ) and the Faint Images of the Radio Sky at Twenty- 
centimetres (FIRST) survey (Bec ker et al.||1995] l provide radio 
measurements at 1.4 GHz (21 cm) for a large fraction of the sky 
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Table 2. Available radio imaging data for some RQQs. 



Name 


v a 


Beam b 


Size c 


PA d 


Reference 


PG 0026+129 


8.4 GHz 


0'.'35 


0'.'60 




(1) 


PG 0052+251 


4.8 GHz 


0'.'59 


< '.'09 




(1) 


Mrk 1014 


8.4 GHz 


0'.'36 


2'.'6 


90° 


(1) 


PG 1612+261 


8.4 GHz 


0'.'28 


l'.'O 


66° 


(2) 


PG 1700+518 


8.4 GHz 


0'.'28 


0'.'5 


-42° 


(2) 


PG 1700+518 


1.6 GHz 


0'.'02 


0'.'5 


-41° 


(3) 


PG 2130+099 


4.8 GHz 


0'.'35 


2'.'9 


-48° 


(1) 



Notes. (a) Observing frequency v of the radio observations. <b) Beam 
size of the interferometric radio observations. (c) Maximum size of the 
extended radio emission. (d) Position angle of the major radio axis. 90° 
points to the East and 0° to the North. 

References. (l) |Leipski et aTJpOOo) ; (2) |KukulaetaTJ(T998) ; (3) |Yangl 
[etaLlpOT2) 



with a completeness level of 2.5 mJy and 1 mJy, respectively. 
We adopt a power-law dependence of the radio flux densities 
(f v oc v a ') with a spectral index of a T = -0.5, at the border line 
between steep- and flat-spectrum radio sources, to estimate the 
flux density at 5 GHz from measurements at 1 .4 GHz when re- 
quired. Only 4 objects in the sample can be classified as RLQs, 
8 objects have intermediate radio properties, and the majority of 
objects within the sample are clearly RQQs. 

Besides the integrated radio measurements, archival high 
spatial resolution images for 6 of our RQQs were obtained with 
the VLA interferometer at subarcsecond resolution and were 



presented by |Kukula et al.| ( |1998| l and |Leipski et al.| (|2006). We 
retrieved the processed images from the VLA image archive^] 
and some have also been kindly provided by the authors. The 
characteristics of these images are listed in Table [2] All except 
PG 0052+251 display extended radio emission with sizes up to 
3". A triple radio source as an unambiguous sign of a jet was 
identified in the RQQs Mrk 1014, in PG 2130+099 using the 
VLA, and in PG 1700+518 using the higher spatial resolution of 
the European VLBI network ( |Yang et al.|2012| l. 

2.3. QSO host morphologies 
2.3.1. Available broad-band imaging 

High spatial resolution broad-band images were obtained with 
the HST for most of the PG QSOs over the last decade, except 
PG 1427+490, PG 1612+261 and PG 2214+139. These HST 
images were initially published in various unrelated articles (see 
Table [3]), but a homogeneous analysis of all the public data was 
presented by Kim et al. ( |2008| l. Ground-based broad-band im- 
ages are additionally available for all QSOs covered by the SDSS 
footprint. The SDSS survey images with an effective exposure 
time of 54s were usually obtained at a seeing around ~1'.'2 (r). 
Yet, it has so far not been attempted to uncover the underlying 
QSO host galaxies in these images given the short exposure time 
and the low spatial resolution of the images. 

No broad-band images have so far been presented for most 
of the HES QSOs in our sample. Ground-based H band images 
were obtained for a few QSOs in 2003 with the SOFI IR spectro- 
graph and imaging camera ( |Moorwood et al.| 1 998) > mounted on 
the New Technology Telescope (NTT) and are publicly available 
in the ESO archive, which we reduced with the corresponding 



https :// archive . nrao . edu/archi ve/ archi veimage .html 



standard ESO pipeline. Additionally, we obtained R band im- 
ages for three HES QSOs with the 60 cm telescope of the Light- 
buckets robotic telescope network (hereafter LB). This small 
telescope is equipped with an Apogee Alta U42 CCD camera 
providing a spatial resolution of 0'.'57 pixeL 1 . With an exposure 
time of ~1 h, these images are nominally deeper than SDSS, but 
suffer from a seeing of >2". We reduced these images following 
the standard procedures including bias subtraction, flat-fielding 
and combination of the dithered exposures with custom Python 
scripts. 

2.3.2. QSO-host deblending 

In order to perform a consistent analysis of the available broad- 
band images we deblended the QSO and host galaxy compo- 
nent using GALFIT ( |Peng et al.l|20T()} [2002] ). We emphasise 
that our main aim here is only to recover the stellar light sur- 
face brightness distribution of the underlying host galaxies and 
distinguish between bulge-dominated and disc-dominated host 
galaxies whenever possible. A more detailed analysis of the im- 
ages would be beyond our scope and hardly possible given the 
low spatial resolution of many images. 

We created an empirical Point Spread Function (PSF) from 
bright stars in the field close to the target. Dedicated PSF stars 
were observed with HST only for a few objects, so that we relied 
on synthetic PSFs generated with the TinyTim software package 
( |Krist|1995] l. First, we used GALFIT to model each object with a 
PSF component only, to check whether any significant residuals 
of a host galaxy could actually be detected. In a second step, we 
modelled the images with a superposition of a Sersic light profile 
( |Sersic|[T968) for the host and a PSF for the QSO component. 
Two examples of the deblending results for the ground-based 
SDSS image are presented in Fig. [2] 

Each QSO host was modelled with three different radial Ser- 
sic surface brightness profiles: 1) Sersic index fixed to n = 1 
(exponential profile), 2) Sersic index fixed to n — 4 (de Vau- 
couleurs profile), or 3) Sersic variable as a free parameter. The 
resulting x 2 values and the best-fit structural parameters are re- 
ported in Table [3] The superposition of a bulge and a disc model 
provides a significantly better fit to the surface brightness dis- 
tribution of PG 2130+099 and I Zw 1 than a single component. 
The structural parameters reported in Table [3]for these two QSO 
hosts refer to the disc component of the model. 

The QSO hosts observed with HST are all well resolved and 
we recover the host galaxies also for the majority of the ground- 
based observations. PSF-subtracted host galaxy images are pre- 
sented below in Fig.|7]for all objects. The ground-based images 
clearly suffer from their low spatial resolution, so that the mod- 
els are strongly affected by PSF mismatches in compact galaxies. 
Such mismatches sometimes lead to non-physically high Sersic 
indices. The comparison of the SDSS i and r images allow us 
to judge the robustness of the results. Strongly different results 
indicate that the structural parameters are not well constrained. 
Furthermore, a close apparent companion affected the modelling 
of SDSS J1655+2146, so that a meaningful model could only be 
inferred for the host galaxy with a fixed Sersic index. 



2.3.3. Morphological classification 

Discriminating between a bulge- and a disc-dominated QSO host 
galaxy is not straightforward, even for the high spatial resolu- 
tion images obtained with HST. A commonly used criterion is 
the best-fitting Sersic index n, where n > 2.5 indicates a bulge- 
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Fig. 2. Two examples of the QSO-host deblending process for SDSS 
i band images. The PSF was constructed from a bright nearby star in 
the field. A best-fitting Sersic surface brightness model plus a PSF for 
the nucleus was obtained with GALFIT and subtracted from the original 
SDSS image to highlight the residuals. The host galaxy is recovered 
after subtracting the best-fit PSF component from the original image 
with nucleus-to-host ratios of -3.5 and -2.2 for SDSS J0057+1446 and 
SDSS J0836+4426, respectively. 



dominated and n < 2.5 indicates a disc-dominated systems (e.g. 
|Bell et al.|2004||Sanchez et al.|2004b| l. Alternatively, the lowest 
X L value is computed for a bulge and a disc-model to identify the 



best matching model (e.g. 


Dunlopetal.2003 Guyon et al. 2006 


Jahnke & Wisotzki|20031 


McLeod & McLeod|2001[|Sanchez &| 


Gonzalez-Serrano 2003). However, the statistical significance of 



the favoured model compared to other options needs to be con- 
sidered. On the other hand, a visual classification of the host 



is often preferred to identify spiral structures of disc-dominated 
systems, or apparently disturbed merging system. 

For the HST images, the x 1 values robustly discriminate be- 
tween the models. We adopt the criterion by Guyon et al. (2006), 
which demands a significant difference in the reduced ;pof 
\xlin = 1) - xl( n = 4)1 > 0.05. The y 2 classification always 
agreed with the one based solely on the best-fit Sersic index and 
we achieved good agreement with literature results. The situa- 
tion is more difficult for the ground-based observations, where 
the^- 2 values sometimes differ only marginally. A visual inspec- 
tion of the images, together with large axis ratios (b/a), may pro- 
vide evidence in some cases for disc-dominated hosts even when 
a bulge-dominated model was formally preferred as in the cases 
of RBS 219, PG 1612+261, PHL 1811 and HE 2353-0420. 
When the best-fitting Sersic index is n > 7, we doubt that any 
reliable classification can be made, because flux transfer likely 
happened between the PSF and the host galaxy component. In 
the other cases, we refer to the best-fitting Sersic indices for our 
classification, but consider these classifications as tentative. 

From the 28 QSOs with available broad-band imaging, we 
classified 17 host galaxies into 5 bulge-dominated, 8 disc- 
dominated and 4 ongoing merging galaxies. 5 of the remain- 
ing hosts are likely bulge-dominated systems and one is proba- 
bly disc-dominated. For the rest of the host galaxies we were 
not able to constrain the morphology. Our sample seems to 
be almost evenly split up into bulge- and disc-dominated sys- 
tems without any preference. A similar balance of morpholog- 
ical types were also reported for completely different samples 
of QSO host galaxies (e.g. Cis ternas et al.||20TTj |Jahnke et al.| 
|2004) . 



3. Spatially resolved spectroscopy of QSOs 

3.1. PMAS IFU observations 

IFU observations for the initial sample were carried out with 
PMAS in the nights from the 2th to 7th of April in 2002 and from 
the 2nd to 4th of May 2003. The QSOs of the extended sample 
were observed in the nights from the 13th to 17th of March and 
from the 5th to 8th of September in 2008. In all cases we used 
the PMAS lens-array with the highest magnification providing a 
seeing-limited sampling of 0'.'5x 0'.'5 with a field of view (FoV) 
of 8" x 8" for the 256 spaxel^] The host galaxies of the ob- 
served QSOs effectively fill the FoV in most cases (see Fig. 17}, 
but are significantly larger for the three nearest objects (z < 0. 10) 
and a few extended systems at higher redshift. 

The low resolution V300 grating was used for the obser- 
vations of the initial QSO sample covering a large wavelength 
range including the major diagnostic spectral features around 
the Hy6 and Ha emission lines simultaneously. Higher spectral 
resolution observations with the V600 grating were taken for the 
extension of the sample, but then centred on the H/3 emission 
line only. For both setups we estimated the spectral resolution 
from the O i A5571 night-sky emission line to be 6. 1 A and 3.2 A 
Full- Width at Half Maximum (FWHM), respectively. At least 
two exposures were usually taken for each QSO with individual 
integration times in the range of 900 s up to 2400 s reaching a to- 
tal integration time of 1-2 hours per object. Table|4]summarises 
the details of the observations including notes on the ambient 
observing conditions. 

2 spectral pixel, a single spatial resolution element containing a spec- 
trum along the entire wavelength range 
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A set of arc lamp (HgNe) and continuum lamp exposures 
were acquired before or after the target frames for fibre trac- 
ing and wavelength calibration. Skyfiats were taken during twi- 
light to measure the variations in the fibre-to-fibre transmission, 
and spectro-photometric standard stars were regularly observed 
at the beginning and at the end of the nights with the same in- 
strumental setups as the QSOs. 

3.2. Data reduction 



more exposures were available to reject remaining cosmic rays 
and artefacts. 

Most of the data were taken under clear or photometric con- 
ditions. In order to check the quality and uncertainty of the 
spectrophotometry, we compared the flux-calibrated QSO spec- 
tra with available SDSS spectra taking into account the size of 
SDSS fibres. The deviation between the QSO spectra in the con- 
tinuum over the whole wavelength range was < 20%, which we 
adopt as our systematic spectrophotometric uncertainty. 



The basic data reduction was performed with the R3D reduc- 
tion package (Sanchez 2006) designed to reduce fibre-fed inte- 
gral field spectroscopy data. We used the independent modular 
tasks of the package for the following basic reduction steps: 1 . 
bias subtraction, 2. fibre tracing on the continuum lamp expo- 
sure, 3. extraction of the fibre spectra, 4. wavelength calibration 
with the arc-lamp exposures, 5. fibreflat correction based on ob- 
served skyfiats, 6. flux calibration of the spectra using spectro- 
photometric standard star observations, and 7. creation of a data- 
cube from the row-stacked spectra. 

Prior to the reduction we detected and removed cosmic rays 
hits from the individual raw frames with the DCR routine (Pych 
|2Q04) > that was designed originally for long-slit spectroscopic 
data. We used a conservative configuration of DCR such that 
night sky lines or emission-line features of our target itself were 
not removed. A few prominent cosmic ray hits inevitably re- 
mained in the data given the imperfection of the algorithm for 
IFU data. Thus, we always required that emission lines must be 
spatially extended, that they have a line width equal or greater 
than the instrumental resolution, and that the doublet line ratios 
are consistent with their theoretical values. 

As discussed in |Rofh et al.] ( |2005[ ), the PMAS instrument 
suffers from flexure depending on the hour angle of the observa- 
tions. This inevitably leads to shifts between different exposures 
of the same target in dispersion and/or cross-dispersion direc- 
tion of the spectra on the charge-coupled device (CCD) detec- 
tor. Flexure shifts were measured with the IRAF task imalign 
by comparing the position of ~20 sky emission lines distributed 
over the whole CCD in the raw image. When only one set of 
calibration frames (continuum+arc lamp) was available for a se- 
quence of observations, we resampled the calibration frames to 
match with position of the individual exposures. The science 
frame taken right after or before the calibration frames was used 
as the reference for the zero position. 

A median sky spectrum was created from spaxels containing 
no emission lines from the host galaxy and subsequently sub- 
tracted from the whole datacube. Given that several QSO hosts 
fill the entire PMAS FoV we accepted a contamination of the sky 
spectrum by stellar continuum emission in these cases. This has 
no strong impact on our emission line analysis, but the Balmer 
emission lines need to be corrected for the underlying stellar ab- 
sorption lines. How the inaccurate sky subtraction affects our 
emission-line measurements will be addressed in Sect. 13.41 

The datacubes were then corrected for differential atmo- 
spheric refraction (DAR) by resampling each monochromatic 
slice such that the position of the QSO is aligned to a com- 
mon reference point at all wavelengths. All standard star ob- 
servations were used to create a mean sensitivity function for 
each run and instrument setup, so that the datacubes are prop- 
erly flux calibrated taking into account the airmass and the at- 
mospheric extinction law at the Calar Alto observatory (Sanchez 



et al. 



20071. Finally, we combined the datacubes using either 



3.3. QSO-host deblending with QDeblend 3D 

Similar to the deblending of the broad-band images we need to 
decompose the QSO and host emission also in our IFU data to 
study and interpret the emission of the host galaxy without the 
QSO contamination. The three-dimensional nature of the IFU 
date requires new analysis techniques and software. We de- 
veloped the software tool QDeblend 3D for this particular tasl^] 
At the core of the program is an iterative deblending algorithm 
which is a significant improvement of a simple technique intro- 
duced by Christensen et al.|([2006 ). 
j i 1 1 i 

The basic concept is to treat the spaxels of the IFU datacube 
as a set of independent spectra rather than as a sequence of 
monochromatic images. In the absence of atmospheric disper- 
sion the spectrum of a point source would then be the same in 
each spaxel, multiplied by a scale factor according to the PSF 
In the case of QSO observations, the scale factors of the PSF 
can be directly determined from the strength of the broad emis- 
sion lines with respect to the adjacent continuum as described 
by pahnke et aL]( |2004| l. Here we outline the iterative algorithm 
implemented in QDeblend 3D . A more detailed explanation with 
illustrative diagrams of the iterative process can be found in the 
user manual of QDeblend 3D . 

In a first step a high S/N nuclear QSO spectrum is extracted 
from the datacube. After a PSF has been determined from a 
broad QSO emission line, the nuclear QSO spectrum is sub- 
tracted from each spaxel after multiplied by the corresponding 
PSF scale factor. However, significant over-subtraction will oc- 
cur when the nuclear QSO spectrum is contaminated by host 
galaxy light. In QDeblend 3D we minimize this host galaxy con- 
tamination by iteratively estimating a mean host galaxy spectrum 
from the residual datacube that is subsequently subtracted from 
the QSO spectrum. A rectangular aperture and surrounding an- 
nulus is used as a default to extract QSO and host galaxy spec- 
tra, but any set of spaxels for the QSO and host galaxy spectrum 
could be specified. Because the host galaxy becomes brighter 
towards the galaxy centre, the mean host spectrum from the an- 
nulus should be scaled in brightness before it is subtracted from 
the initial nuclear QSO spectrum. That scale factor is uncon- 
strained from the IFU data itself. One may either adopt an a- 
priori defined analytic surface brightness profile, a constant sur- 
face brightness, or a manually set scale factor which fulfils a cer- 
tain criterion. The iteration is stopped when the decontaminated 
nuclear QSO spectrum converges. 

For our PMAS spectra we use the broad Ha or H/3 lines of 
the QSOs to perform a QSO-host deblending with QDeblend 3D . 
Because we focus on the emission lines of the ionised gas 
around Hfi and Ha region in this work, we restrict the QSO-host 
deblending to this narrow wavelength regions around the two 
Balmer lines. This almost entirely avoids the problem of PSF 
mismatch due to the wavelength dependence of the atmospheric 
seeing. The nuclear QSO spectra are taken from the brightest 



a variance weighted mean or a clipping method whenever 3 or 3 available for download at http://sf.net/projects/qdeblend/ 
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Fig. 3. Comparing the PMAS spectra of PKS 2349-014 (region A) 
with two SSP model spectra matched to the HST photometry. The sys- 
tematic offset between the SSP models and the PMAS spectra is caused 
by an over-subtraction of the sky background in the PMAS observa- 
tions. 
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Fig. 4. Underestimation of the H/3 emission line flux (A log(H/3)) due 
to the stellar absorption line as a function of the stellar population age 
for the PMAS spectrum of PKS 2349-014 (region A). 



spaxel in the data, and the surrounding 8 spaxels are used to es- 
timate a mean host galaxy spectrum from the residuals. Here 
we assume a constant surface brightness to extrapolate the host 
galaxy emission towards the central spaxel. The QSO spectra 
converged already after four iterations at which we stopped the 
iterative process. 



teresting to study the stellar populations of the QSO hosts from 
their continuum emission, we find that given our exposure times, 
the S/N is generally too low to quantitatively interpret the stellar 
absorption patterns. Nevertheless, the measured strength of the 
Balmer emission lines will depend on the depth of the underly- 
ing stellar absorption, which we take into account when using 
the Balmer lines for diagnostic purposes. In order to explore the 
impact of the Balmer absorption line on our measurements, we 
combined synthetic single stellar population (SSP) models from 
Bruzual & Chariot ( |2003| l with the absolute surface brightness 
photometry of the HST broad-band images. In a first step, we 
convolved the nucleus-subtracted HST images with the PSF of 
the corresponding PMAS observation and rebinned the image to 
the spatial sampling of our PMAS data. Afterwards, we scaled a 
suite of SSP spectra to the surface brightness in several specific 
regions. In Fig. [3] we compare two SSP models (300 Myr and 
5Gyr) matched to the HST photometry with the PMAS spec- 
trum of PKS 2349-014 (region A as defined in Fig. [7]i as an 
example. The PMAS spectrum is well below both SSP model 
spectra, which can be entirely attributed to the over-subtracted 
sky background. 

The strengths of the Balmer absorption lines depend on the 
age of the stellar population. We measured the differences be- 
tween the H/3 emission line fluxes of the continuum subtracted 
and unsubtracted PMAS spectra expressed as Alog(H/3) = 
log(Hj6 su b.) - log(Hy6 unsu b.) for 18 SSP spectra with ages of the 
order of 300 Myr < fssp < 5 Gyr assuming solar metallicity. We 
show the dependence of AH/3 on f SSP for PKS 2349-014 (region 
A) as an example in Fig. [4] We find that the H/3 emission line 
flux is underestimated due to the corresponding stellar absorp- 
tion line by 0.07-0.2 dex in this case of PKS 2349-014. The 
redshift of that object (z = 0.174) is close to the median redshift 
of our sample and representative for the majority of objects. The 
H/3 emission-line flux is more severely affected for the few object 
at z = 0.06 given the lower equivalent width of their emission 
line. 

Because of the unknown age of the stellar population, we 
used the average of Alog(HyS) taken over all SSP ages with an 
error corresponding to the possible range of values. For the 
host that fill the PMAS FoV and were not observed with HST, 
we corrected all H/3 lines by an average factor of A log(H/3) = 
0.13 ± 0.08 dex, determined from all objects with reliable HST 
photometry. Although the Ha emission line is also affected by 
the stellar absorption, we made no correction given that the cor- 
rection factor (0.03 + 0.01 dex) is negligible for our purposes. 



3.4. Imperfect sky subtraction and Balmer absorption 

One problem with our PMAS IFU observations is the imper- 
fect subtraction of the sky background. Because no blank sky 
fields were observed together with the target, we were forced to 
estimate the sky directly from the target datacubes itself. We 
used spatial regions free of extended emission lines to avoid any 
over-subtraction of these lines, but those regions may contain 
significant emission from stellar continuum if the host galaxy 
extends across the PMAS FoV. Thus, the stellar continuum is 
over-subtracted in several cases for which we were unable to re- 
liably infer the "real" stellar continuum from our data. 

Our PMAS observations were primarily targeted to map the 
emission lines of the ionised gas. While it would be also in- 



4. AGN parameters from the QSO spectra 

Our PMAS data provide high S/N spectra of the QSOs without 
any wavelength-dependent slit losses, because the IFU covers 
the entire seeing disc of a point source even under bad seeing 
conditions. We use those spectra to derive some fundamental 
AGN parameters, like redshifts, black hole masses Mbh, and 
Eddington ratios A = Lboi/^Edd- The contamination of the QSO 
spectrum by host galaxy light was removed or at least minimised 
by our QSO-host deblending process described above. For a few 
QSOs where the SDSS spectra had significantly higher S/N we 
used these instead. From the SDSS broad-band images we es- 
timated a host contamination of <15% in the SDSS spectra and 
corrected the continuum luminosities correspondingly. 
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rest-frame wavlength [A] rest-frame wavlength [A] 

Fig. 5. Left panel: QSO spectrum of PHL 1811 observed with PMAS as a representative example for a strong Fen emitter. The available 
wavelength range of the V300 grating was set such that H/? and Ha would be covered at the redshift of the QSO. The Fe n complex blueward of 
H/3 is not fully covered for PHL 1811, as in many other objects. Right panel: Comparison of two available Fen templates from Boros on"& Green| 
(1992) and Veron-Cetty et al. (2004) both based on the proto-typical QSO I Zw 1. Both templates were scaled to match the isolated Fen /15014 
line. The disagreement between the templates is severe given that they were generated from the same QSO. 



4.1. Spectral modelling of the H/3-[0 \\\] region 

All of the AGN parameters of interest were calibrated to spec- 
tral features around the broad H/3 line. This wavelength region 
is bracketed by two adjacent Fen complexes as shown in Fig. [5 
(left panel). Most studies o f QSO spectra (e.g.lBoroson & G reen 
[T992l|Greene & Ho|2005[|McGill et al.|2008[ |Woo et al.|2006] > 
used a template generated from the proto-typical QSO I Zw 1 
( |Boroson &Green|1992t|Veron-Cetty etal.|2004| l to subtract the 
Fe n complexes over a large wavelength range. We cannot apply 
this Fe ii template fitting method here in a consistent manner be- 
cause most of our QSO spectra do not reach enough to the blue 
to reliably estimate the underlying power-law continuum. We 
also note that the two available Fe n templates significantly dis- 
agree (Fig. [5] left panel). However, the dominant Fen emission 
in the proximity of H/3 and [O m] originates in both cases from 
two isolated Fen lines, enabling us to choose a more simpler 
approach to handle the Fe n /Li4929, 5018 lines over a restricted 
wavelength range. 

We thus modelled only the wavelength region of 4700- 
5100A, which encompasses H/3, [Oin] and almost pure contin- 
uum regions on either side. A superposition of Gaussians was 
used to represent the line profile of the broad and narrow emis- 
sion lines including the two prominent Fe II/LI4929, 5018 lines. 
The full model consists of up to 3 Gaussian components for the 
H/3 line, up to 2 Gaussian components for the profile of the [O in] 
and Fen doublets and a straight line to approximate the local 
continuum. The two lines of a doublet are coupled as a sys- 
tem in line width, radial velocity and line ratio (a factor of 3 for 
[Om] X14960,5007 and 1.29 for Fen M4924, 5018), which sig- 
nificantly increases the robustness of our model and reduces the 
number of free parameters. If a narrow H/3 component is present 
on top of a broad H/3 emission profile, we couple its kinematic 
parameters with that of the narrow [Om] component as these 
lines originate both from the NLR. In cases where the Fe n lines 
were broad and their width could hardly be constrained against 
the background continuum, we coupled the line width and radial 
velocity to that of the broad Hp line. 

The best-fitting parameters were found by searching for the 
lowest x 2 of the residuals with a downhill simplex algorithm 
(e.g. |Press et aT7|| 1992) >. Two representative examples of these 
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Fig. 6. Spectral modelling of the H/?-[0 in] region for two representa- 
tive QSO spectra. The black line correspond to the observed QSO spec- 
tra and the red curve represents our best fit model. The model consists 
of individual Gaussian components for the emission lines (H/3-magenta, 
[O m]-green, Fe n-blue) above a straight line for the local continuum. 
The residual of the model are shown in a dedicated panel below. 

models are presented in Fig. [6] From the best-fitting model, 
we measured several quantities including the continuum flux at 
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Table 5. AGN parameters measured or derived from the QSO spectrum 



Object /NLR,[Om] l0gi N LR,[0m] /5100 logi 5 100 FWHMtyj 0-^8 log(M BH ) log A 

[ lr ' 6 3] [ er g s_ '] [ lr ' 6 ^l] [ er g' s_1 ] [ kms_1 ] [ kmsr '] [ M °J 



PG 0026+129 


505 


±3 


42.54 


±0.09 


27.5 + 


6.7 


44.96 ± 


0.11 


2430 ± 


28 


3471 


± 18 


9.0 


-1.2 


IZw 1 


1667 


± 13 


42.23 


±0.09 


45.8 ± 


11.3 


44.31 ± 


0.11 


1111 ± 


12 


1652 


± 10 


8.0 


-0.9 


PG 0052+251 


667 


±2 


42.69 


±0.09 


25.1 + 


5.8 


44.98 ± 


0.10 


4604 ± 


24 


2088 


±8 


8.6 


-0.7 


SDSS J0057+1446 


203: 


t 14 


42.27 


±0.12 


18.6 + 


4.2 


44.95 ± 


0.10 


10506 ± 


159 


4143 


±39 


9.1 


-1.3 


HE 0132-0441 


200: 


t33 


42.16 


±0.16 


12.6 + 


2.9 


44.69 ± 


0.10 


1487 ± 


21 


960: 


t47 


7.7 


-0.2 


HE 0157-0406 


82: 


t2 


42.08 


±0.10 


6.2 ± 


1.3 


44.72 ± 


0.09 


5493 ± 


174 


2333 


±77 


8.5 


-1.0 


SDSS J0155-0857 


211 : 


t 11 


42.23 


±0.11 


6.7 ± 


1.4 


44.47 ± 


0.09 


6650 ± 


87 


2824 


±35 


8.6 


-1.2 


Mrk 1014 


901 : 


t 10 


42.85 


±0.09 


13.7 + 


3.0 


44.77 ± 


0.09 


3002 ± 


25 


1549 


± 16 


8.2 


-0.6 


SDSS J0846+4426 


1038 


±53 


43.34 


±0.11 


17.9 + 


3.9 


45.35 ± 


0.10 


4169 ± 


128 


1771 


±56 


8.6 


-0.4 


SDSS J0948+4335 


179 


±5 


42.45 


±0.10 


5.4 ± 


1.2 


44.72 ± 


0.09 


4060 ± 


61 


1571 


±30 


8.2 


-0.6 


SDSS Jl 131+2623 


252: 


t33 


42.68 


±0.14 


3.4 ± 


0.7 


44.59 ± 


0.09 


9370 ± 


328 


3957: 


t 128 


8.9 


-1.5 


PKS 1217+023 


281 


±2 


42.71 


±0.09 


14.6 + 


3.2 


45.20 ± 


0.10 


4300 ± 


178 


2272: 


t686 


8.7 


-0.7 


3C 273 


996: 


t32 


42.87 


±0.10 


272.7 + 


59.2 


46.04 ± 


0.09 


2852: 


t7 


1999 


±7 


9.1 


-0.2 


SDSS J 1230+ 11 00 


241 


±7 


42.64 


±0.10 


5.1 ± 


1.1 


44.73 ± 


0.10 


4439 ± 


129 


2307 


±30 


8.5 


-0.9 


HE 1228+0131 


454: 


t27 


42.23 


±0.11 


44.5 + 


9.7 


44.95 ± 


0.09 


1671 ± 


11 


1572 


±38 


8.3 


-0.5 


SDSS J1230+6621 


447 


±7 


42.65 


±0.09 


5.9 ± 


1.2 


44.54 ± 


0.09 


3297 ± 


93 


2190 


±22 


8.4 


-1.0 


PG 1427+480 


240 


±5 


42.55 


±0.10 


5.3 ± 


1.1 


44.67 ± 


0.09 


2455 ± 


68 
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8.0 


-0.4 
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-1.0 


HE 1453-0303 
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42.53 


±0.12 


27.4 + 


7.9 


45.31 ± 


0.12 


2513 ± 


43 


1795 


±27 


8.6 


-0.4 


PKS 1545+21 


286 


±2 


42.84 


±0.09 


14.1 + 


3.1 


45.29 ± 


0.10 


5644 ± 


38 


2397 


± 15 


8.8 


-0.7 


P(j 1612+261 


2591 


+ 17 


43.14 


± 0.09 


21.2 + 


5.1 


44.74 ± 


0.1 1 


2770 ± 


26 


1 177 


± 11 


7.9 


-0.4 


SDSS J1655+2146 


687 


±7 


42.71 


±0.09 


12.0 + 


3.0 


44.65 ± 


0.11 


2483 ± 


37 


1576 


±27 


8.1 


-0.6 


PG 1700+518 


145 


±5 


42.62 


±0.10 


15.6 + 


3.4 


45.43 ± 


0.09 


2067 ± 


41 


2600 


±53 


9.0 


-0.7 


PG 2130+099 


870: 


t 15 


41.98 


±0.09 


48.3 ± 


11.1 


44.40 ± 


0.10 


2727 ± 


10 


2144 


± 17 


8.3 


-1.0 


HE 2152-0936 


122 


±7 


42.17 


±0.11 


60.7 ± 


13.8 


45.59 ± 


0.10 


1877 ± 


13 


1853 


± 16 


8.8 


-0.3 
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241 


± 1 


42.59 


±0.09 


9.6 ± 


2.3 


44.89 ± 


0.11 


5486 ± 


45 


2325 


± 19 


8.6 


-0.9 


HE 2158+0115 


124 


± 1 


41.92 


±0.09 


4.8 ± 


1.1 


44.21 ± 


0.10 


2739 ± 


33 


1637 


±42 


7.9 


-0.9 


PG 2214+139 


466: 


t 11 


41.78 


±0.10 


48.3 ± 


12.1 


44.44 ± 


0.11 


5506 ± 


23 


3190 


±77 


8.6 


-1.4 


HE 2307-0254 


61 : 


t3 


42.00 
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8.6 ± 


1.9 


44.89 ± 


0.10 


2172 ± 


27 


1823 


±32 


8.4 


-0.7 


PKS 2349-014 


279 


±7 


42.40 


±0.10 


14.0 + 


3.0 


44.85 ± 


0.09 


4041 + 


83 
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±25 


8.6 


-0.9 


HE 2353-0420 


224 


± 1 


42.58 


±0.09 


2.7 ± 


0.6 


44.42 ± 


0.10 


4176 ± 


40 


1774 


± 17 


8.1 


-0.9 



5 100 A (/5100) and the width of the broad H/3 component ex- 
pressed both in terms of FWHM (FWHMh/j) and line dispersion 
(°~Hfi)- These values are listed in Table [5] for the entire sample. 
To check whether these measurements are robust compared to a 
full Fe n template fitting approach we performed a simple con- 
sistency check. We generated artificial QSO spectra from our 
best-fitting models excluding the Fen components but adding 
the |Boroson & Green| iron template matched to the Feii/15014 
line strength and broadened to the FWHM of the broad H/? line. 
These artificial spectra were subsequently analysed in exactly 
the same way as the original QSO spectra. The recovered values 
for the FWHM of H/3 and the continuum flux at 5 100 A differ by 
only 0.02 dex in the mean, with a maximum offset of 0.05 dex 
and 0.08 dex, respectively. 



4.2. Black hole masses and Eddington ratios 

Black hole masses can be estimated from single epoch broad- 



line AGN spectra using the so-called virial method ( Peterson & 
|Wandel||2000| Vestergaard 2002). It combines the empirically 
derived size-luminosity relation of the BLR as determined via 
reverberatio n mapping (e.g.|Bentz et al.|200 9 2006 ; Kaspi et al. 
[20051 |2000l |Peterson et al.|2004| > and the kinematics of the BLR 
estimated from the width of the broad emission-lines assuming 
virial motion of the clouds. We adopted the empirical relation 
between the BLR size and the continuum luminosity at 5100 A 
(£5100) calibrated by Bentz et al. ( 2009} , and the prescription 



of |Collin etaT] ( |2"006 



to infer the kinematics of the broad line 



clouds from the H/3 line dispersion (0- ^3). A virial scale factor 
of 3.85 was estimated by |Collin et al. (2006! to match the virial 
BH masses with the Mbh - cr* relation of local inactive galaxies 
( |Tremaine et al.|2002| l. This yields a black hole mass calibration 
of 



Mr 



10' 



( 



0"H/3 



L 5 



100 



\ 1000 km s- 1 / VlO^ergs- 1 



0.52 



(1) 



The error on the estimated BH masses using the virial method 
is dominated by systematic errors rather than measurement er- 
rors from high S/N QSO spectra (see |Denney et al.] 2009| >. In 
addition to the intrinsic uncertainties in the adopted scaling rela- 
tions and the adopted line profile model, the unknown geometry 
and kinematics of the BLR and the correct determination of the 
QSO continuum flux are sources of uncertainty. The BH masses 
estimated in this way (see Table [5]) are assumed to be accurate 
within a factor of ~2 (~0.3 dex). 

Furthermore, we computed the bolometric luminosity of the 
QSOs adopting a standard bolometric correction of 



9x^(5100A) 



(2) 



following Kaspi et al. (2000) which is in agreement with the es- 
timates of Richards et al. (2006). From the BH masses we also 
calculated the Eddington luminosities Lsdd and Eddington ratios 
A = Lboi/^Edd- The derived values for the Mbh and A are given 
in Table [5] Almost all of our QSOs have black hole masses of 
Mbh <: 1O 8 M and Eddington ratios in between 0.1 < A < 1. 
None of our QSOs accretes above the Eddington limit, which is 
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in agreement with the results of low-redshift QSO studies (e.g. 
ISchulze & Wisotzki|[20T0l|Steinhardt & ETyisPJflty . The me- 
dian BH mass and median Eddington ratio for our sample are 
M B h - 4 x 10 8 M G and A - 0.15. In individual cases, M B h are 
found to be 0.5 dex larger than previous estimates based on the 
FWHM of H/3 (e.g. |Vestergaard & Peterson|2006| l , but the values 
based on cr ap should be more robust as argued by |Peterson et al. 
pOQ4l > and |Collin et aE1 ( |2006| >. 

5. Emission line diagnostic of extended ionised gas 

5. 1 . Mapping of extended emission and their spectra 

After subtracting the QSO contribution with QDeblend 3D from 
the observed data we study the properties of ionised gas around 
the QSOs. Firstly, we extract 40 A wide narrow-band images 
centred on the [Om] and Ho- lines of the QSO after remov- 
ing any continuum emission estimated the adjacent bands. The 
synthetic passband of our narrow-band images roughly corre- 
spond to ±900kms~ 1 in the rest frame at the median redshift 
of the sample, so that extended emission even with high ve- 
locities would be captured. We detect extended emission in 
[Om] around 16 of 31 QSOs (-50%), of which a high fraction 
have also detectable Ha emission (when covered by our obser- 
vations). Whether the apparent [Om] emission around 3C 273 
is a real feature is currently unclear as discussed below. Addi- 
tionally, Ha emission is detected around three objects without 
noticeable [Om] emission, I Zw 1, RBS 219 and PHL 1811. 

In Fig.|7]we present the PMAS narrow-band images of all ob- 
jects for which we detected extended ionised gas in comparison 
to the nucleus-subtracted broad-band host images, when avail- 
able. The ionised gas generally does not extend beyond the size 
of the host galaxies at the surface brightness limits of our ob- 
servations. The only exception is HE 2158-0107 for which the 
ionised gas region is much more extended than the host galaxy 
and reaches even beyond the PMAS FoV with a projected size of 
>22kpc. A detailed discussion and interpretation of this inter- 
esting QSO is presented in Husemann et al. (201 1 1 and beyond 
the scope of this paper. 

As the H/3, Ha, and [N n] lines are rather weak in individual 
spaxels we decided to co-add several spaxels to increase their 
S/N. Based on the morphology of the extended emission we de- 
fined up to 4 characteristic spatial regions, avoiding spaxels with 
obvious signs of artefacts or QSO residuals that could affect the 
emission line measurements. The region boundaries are shown 
in Fig.|7]together with the corresponding co-added spectra in the 
rest-frame wavelength intervals 4800-5 100A and 6450-6750A, 
if covered by our observation. The H/3 and [Om] emission lines 
were clearly detected in most of these objects. Essentially no 
contamination of the broad H/3 emission was left in these ex- 
tended spectra, which demonstrates the robustness of the de- 
blending process at least down to 0'.'5 distance from the QSO 
position. The Ha and [N n] lines are often redshifted into one of 
the atmospheric absorption bands, which prevent a detection of 
these lines in those cases. 

To quantitatively measure the observed emission line prop- 
erties, we modelled each line in the spectra with a Gaussian pro- 
file taking into account the theoretically fixed line ratios of the 
[O in] X14960, 5007 and [Nn] ,Lt6548,6583 doublets. Further- 
more, we assumed that the lines are kinematically coupled hav- 
ing the same redshift and line width. This may not be true in all 
cases, because distinct ionised gas clouds with different kinemat- 
ics and ionisation states can be aligned along our line-of-sight. 
The seeing may enhance this apparent confusion and in a few 



Table 6. Surface brightness limits (3<x) for HST [O m] narrow-band im- 
ages. 



Name 


HST 


smoothed 11 HST PMAS 




[io- 


15 -1 -2 -2l 

erg s cm - arcsec _ 


PG 0026+129 


10.73 


3.79 0.14 


PG 0052+251 


2.77 


0.80 0.18 


Mrk 1014 


2.74 


0.91 0.21 



Notes. (a) A median filter of 5 x 5 pixel box size was used. 



cases the redshift of the extended Balmer and [Om] lines do not 
agree in some regions of the host galaxy. Measured line fluxes, 
line widths and the rest-frame velocity offset with respect to the 
redshift of the narrow QSO [O m] emission line are reported in 
Table [7] Statistical errors on the flux measurements and 3<x up- 
per limits for undetected emission lines were obtained based on 
the noise in the adjacent continuum. 

On two EELR objects we have to comment in more detail: 
The H/3 line flux for the region A and B of PG 1612+261 is un- 
reliable due to a bad column on the PMAS CCD, but regions C 
and D remain unaffected. Off-nuclear long-slit spectra presented 
already by |Boroson et al.| ( |1982| l showed that the [O m]/H/3 line 
ratio is ~10 out to several arcseconds away from the nucleus. 
Thus, we assume that the gas is at least as highly ionised closer 
to the nucleus, although we were not able to reliably constrain 
the line ratios in that region. Furthermore, we find that the [O m] 
emission line in region A is significantly asymmetric requiring 
two independent Gaussian components to model the line profile. 
The QSO 3C 273 actually appears to be surrounded by some 
ionised gas, but it is by far the brightest object in our sample and 
was observed under poor seeing conditions. To check whether 
this feature may simply be an artefact of our data reduction or 
analysis we analysed the 3 individual exposures obtained for this 
target separately. The apparent structure could only be recovered 
in 1 of the 3 exposures implying that this may not be a real detec- 
tion. Until observation with higher quality exists for 3C 273 we 
will treat it as a non-detection throughout the rest of this article. 

No extended line emission was detected around the remain- 
ing 1 1 objects in the sample. We present their narrow-band im- 
ages and available broad-band images in Fig.[8]for completeness. 

5.2. Comparison with HST narrow-band imaging 

Three QSOs among our sample are within the sample of the HST 
narrow-band snapshot survey by |Bennert et al.| (|2002). This en- 
ables us to directly compare the narrow-band images obtained 
from both telescope (Fig. [9). We applied a median smoothing 
filter with a box size of 5x5 pixels to the HST images in or- 
der to catch lower surface brightness features. The 3<x surface 
brightness limits for the HST and PMAS narrow-band images 
are listed in Table[6] The estimated surface brightness limits sup- 
port our visual impression from the comparison of the narrow- 
band images that the PMAS images are at least a factor of 4 
deeper than the HST images even after smoothing. 

PG 0026+129 was observed with the Planetary Camera (PC) 
that has a sampling of O'.'OSpixeF 1 instead of O'.'lpixeF 1 com- 
pared to Wide Field Channel (WFC). Thus, the PMAS image 
for this particular object is roughly 20 times deeper than the 
one of HST, which was most likely the reason why the extended 
[O m] emission is completely missed by HST. The lowest surface 



Article number, page 10 of 



34 



B. Husemann et al.: The properties of the extended warm ionised gas around low-redshift QSOs 




6 4 2 -2 -4 -6 

A« [arcsec] 



- 3 



1 r 1 1 1 1 1 " ■ i 1 1 " i 1 



HST WFPC2 (FU11W) 



























HST WFPC2 (F606W) 
i i i i I ' 



A« [arcsec] 



■ ■Ml 1 II 


| IlltJII 1 1 | llll|M 1 1 














; LJ3001 (« 

' — 1 .... 1 — 1 — 1 




5 10 


5 -5 


-10 -] 




A« [arcsec] 




.... 1 


i i i | i i i 










SDSS M 
' . . . . 1 . . . . 1 . . . 


i .... ' 


5 





-5 -] 




A« [arcsec] 





3 2 1 0-1-2 -3-1 
Aa [arcsec] 



4850 4900 4950 5000 5050 
rest-frame wavelength [A] 



6500 6550 6600 6650 

rest-frame wavelength [A] 





15 10 5 -5 -10 -15 
A« [arcsec] 



3 2 1 0-1-2 -3-4 
Aa [arcsec] 



4850 490O 4950 5000 5O50 
rest-frame wavelength [A] 



6500 6550 6600 6650 

rest-frame wavelength [A] 



I ' I 1 I 1 I ' I 

5kpc 




3 2 1 0-1-2-3 -4 
Aa [arcsec] 



4850 4900 4950 5000 5O50 
rest -frame wavelength [A] 



6500 6550 6600 6650 

rest-frame wavelength [A] 




3 2 1 0-1-2-3-4 
Aa [arcsec] 



4850 490O 4950 5000 5O50 
rest-frame wavelength [A] 



6500 6550 6600 6650 

rest-frame wavelength [A] 




1 


i i i | ■ ■ i , | 


i 1 " 


: .... 1 


... i .... i 


B 

... 1 : ... 1 ... . 



Aa [arcsec] 



4850 4900 4950 5000 5050 
rest-frame wavelength [A] 



Fig. 7. Overview of the individual objects with extended emission after the point-like emission of the QSO was subtracted. Left panels: Nucleus- 
subtracted broad-band images with overplotted contours of extended line emission within the PMAS FoV (black box). A scale bar represents 5 kpc 
at the object redshift and the image is oriented such that North is up and East is to the left. Middle panels: Nucleus-subtracted [O m] narrow-band 
PMAS images with a logarithmic scaling in units of 1CT 16 erg s~'cmr 2 arcsec~ 2 . Note that in the case of I Zw 1, RBS 219 and PHL 1811, the Ha 
narrow-band images are shown instead. The blue boundaries indicate the spatial regions defined to obtain the co-added spectra shown in the right 
panels. A white cross marks the position of the QSO and the black ellipse indicates the FWHM of the PSF. Right panels: Co-added spectra of the 
specific spatial regions around the H/? and Ha line, if it was covered by the spectral range of the used instrument setup. The different spectra are 
arbitrarily scaled and offset from each other, for which the dashed lines indicate their corresponding zero flux density. Vertical dotted lines mark 
the position of emission lines at the QSO redshift, estimated from the peak of the narrow [O in] line or from the peak of the H/3 line when the 
[O in] line was too weak or exceptionally broad. The arbitrarily scaled QSO spectrum is shown as a red line for comparison. Telluric absorption 
line features are highlighted by a crossed circle symbol. 
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Fig. 7. Continued. 
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Fig. 7. Continued. 
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brightness filaments around the nucleus in the HST image may 
tentatively correspond to the high surface brightness structures 
in our PMAS images. It is obvious that |Bennert et al. (2002 ) 
underestimated the size of the extended [O ni] emission for this 
object at least by factor of 3. 

The HST images of PG 0052+251 and Mrk 1014 obtained 
with the WFC are deeper so that high surface brightness struc- 
tures of extended [O m] emission could be detected. Neverthe- 
less, our ground-bases PMAS narrow-band images reach a lot 
deeper, in particular in the case of Mrk 1014 where diffuse emis- 
sion is present between the high surface brightness structures. 
Similar IFU observation of Mrk 1014 obtained by Fu & Stockton 
(2008) display the same [Om] light distribution as our PMAS 
observation suggesting that the structures are real. The low sur- 
face brightness regions south-east of PG 0052+251 was missed 
by |Bennert et al.| ([2002 ) and became visible only after smooth- 
ing the HST data. The maximum [O in] size was thus underes- 
timated by a factor of ~2 in these two cases. We return to the 
issue of NLR sizes in section [6] 

For PG 0052+215 we note further that the contours of the 
smoothed HST images do not agree with the [O m] light distri- 
bution in our PMAS images north of the nucleus. We carefully 
checked that this is not introduced by our data reduction or data 
analysis. A possible origin for the rather symmetric diffuse light 
in the HST image around the nucleus might be continuum emis- 
sion of the host. Bennert et al.| (2002 ) subtracted only the con- 
tribution of the QSO using a PSF-star image, because the S/N 
in the much shallower continuum image of the host was insuf- 
ficient. The unavailability of a continuum image also hampered 
the analysis of the narrow-band HST image of PKS 1545+210 



(3C 323.1) as part of large HST snapshot program to study ex- 
tended gas around radio galaxies ( Privon et al.|2008 ). 

We conclude that HST is of course much better in terms of 
spatial resolution than any optical ground-based observations, 
but the narrow-band images strongly suffer from its low sensi- 
tivity which has already been noted by several authors (e.g. |Ben^| 
|nert|2005"||Privon et al. 2008 ). Additionally, narrow-band images 
require off-band images which are often not sufficiently deep or 
absent to subtract the stellar continuum contribution from the 
on-band images. 

5.3. lonisation sources of the extended ionised gas 

Narrow-band imaging observations of large AGN samples have 
been the main tool for detecting and characterising EELRs (e.g 
Baum et al.||19"8"8"l jBennert et al.||2002| |Mulchaey et al.||1996a[ 



Privo n et al.||2008l |Schmitt et al.||2003a|bl |Stockton & MacK-| 
|enty|1 987). Most of these projects focused on a single emission 
line, which does not allow to constrain the origin of the gas ion- 
isation. Three different ionisation mechanisms for the ISM on 
galactic scales need to be distinguished: 1. Photoionisation by 
the hard UV radiation field. 2. Shock ionisation due to winds 
or radio jets. 3. Photoionisation due to the ionising UV pho- 
tons of young and massive OB stars in star-forming region, i.e. 
classical H n regions. It is possible to distinguish between differ- 
ent excitation mechanisms with the aid of diagnostic emission- 
line ratio diagrams (e.g. Alle n et al.|19 98; Bald win et al.| 1981 ; 
IVeilleux & Osterbrock||1987[ |Vill ar-Marti n et al.||1997| )' The 
most commonly used diagnostic diagram in the optical is the 
[Om] A5007/H/3 vs. [Nn] ^6583/Ha diagram ( |Veilleux & OT] 
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Fig. 8. Overview of the individual objects without detected extended emission. For each object the nucleus-subtracted broad-band images is 
presented, if available, in the left panel with the PMAS FoV indicated by the black rectangle. The black scale bar indicates a physical size of 5 kpc 
at the redshift of the object. The nucleus-subtracted [O in] narrow-band images is presented in the right panels with the a white cross indicating 
the QSO position and the black ellipse indicates again the FWHM of the PSF 
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□ 

Fig. 9. Direct comparison of PMAS and HST [O m] narrow-band images. PMAS narrow-band images of PG 0026+ 1 29 (left panel), PG 0052+25 1 
(middle panel) and Mrk 1014 (right panel) are shown in units of 10~ 16 erg s" 1 cnT 2 arcsec~ 2 . The corresponding HST narrow-band images of 
Ben nert et al.| j2002 1 are overplotted as blue contours also in a logarithmic scaling. 
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Fig. 10. The standard BPT [Om]^5007/H/3 vs. [Nii]^6583/Hq- 
diagram. The 2D histogram of the line ratios for more than 
800 000 emission-line galaxies taken from the SDSS DR7 value-added 
MPA/JHU galaxy catalogue are shown in the grey (logarithmic) scale 
image. All EELR regions of our sample for which both line ratios could 
be measured, even if one ratio is just a limit, are plotted as black sym- 
bols. Demarcation curves by Kewley et al. (2001), Kauffmann et al.| 
( (20031 and[Stasihska eTaTH2008) are drawn as solid, dashed and dotted 
blue lines, respectively. 



terbrock|T9 87), also known as the BPT diagram (Bal dwin et al.| 



1981 1. Those lines are usually the strongest in the optical spec- 
trum and their ratios are most insensitive to reddening. Unfor- 
tunately, AGN photoionisation and shock-ionisation lead to very 
similar line ratios in the BPT diagram and cannot be reliable dis- 
tinguished just from those lines. Anyway, both of these ionisa- 
tion mechanisms are likely related to the AGN and we consider 
all corresponding line ratios being AGN ionised. AGN ionised 
and H n regions populate two distinct branches on this BPT dia- 
gram as shown in Fig. [10] A theoretical limit for the emission- 
line ratios of H n regions that can be reached during a maximum 
starburst event was derived by |Kewley et al.| ( |2001| l as 



Based on a large sample of SDSS emission-line galaxies, |Kauff-| 
|mann et ah] ( |2003| l found that the AGN branch in the BPT dia- 
gram separates from H n branch well below the Kewley et al. de- 
marcation curve and empirically defined a new line that tightly 
encloses the H n branch being parametrised as 



log([Oin]/Hj8) = 0.61/(log([Nn]/Ha) - 0.05) + 1.3 



(4) 



log([Om]/H/3) = 0.61/(log([Nii]/Ho') - 0.47) + 1.19 



(3) 



The region in between the two demarcation lines is often referred 
to as a composite region (Kauffm ann et al.|2003||Kewley et al.| 
2006 ) where the emission-line gas could be partially ionised by 
both mechanisms. However, this classification approach may be 
oversimplified as pointed out by |Cid Fernandes et al.| ( |2010| , be- 
cause even emission-line ratios in the AGN area above the Kew- 
ley et al. curve may have some hidden contribution from Hn 
region. Conversely, a region below that curve could be entirely 
powered by a massive starburst without requiring an AGN con- 
tribution. Furthermore, the AGN branch is composed of Seyfert 
type emission lines and Low Ionisation Nuclear Emission Re- 
gions (LINERs, Heckm an|1980|l with lower [O m]/Hy6 at a given 
[Nii]/Hck ratio. Stasiris ka et al.| ( |2008| ) defined a demarcation 
line between the two classes at a angle of 59°, counter clock- 
wise from a line of constant [O m]/H/J, going through the point 
(log [Nil] /Ha = -0.43, log [Ora]/H/3 = -0.49). As can be seen 
from Fig. [10] LINER-type line ratios are not found for any ex- 
tended ionised gas region in our sample. 

We created a BPT diagram (Fig. [T0| for all the co-added 
spectra in which we could measure both line ratios, [Om] 
A5007/H/3 and [Nn] /16583/Ha. In the case that either H/3 or 
[Nil] /16538 was below the detection limit we used the 3cr up- 
per limit to constrain at least an upper/lower limit on the corre- 
sponding line ratio. Most of the emission-line ratios are on the 
AGN side above the Kewley et al. line for which we assume that 
AGN ionisation is the dominant process to ionise the gas. Re- 
gion C of PG 0052+251, region A of PG 2130+099 and region 
A of PHL 1811 are almost directly on the Hn branch and rep- 
resent regions of ongoing star formation. In case of PHL 1811 
the star forming region is coincident with spiral-like structure 
detected in the high-resolution HST image that was reported by 
Jenkins et aL] ( [20U5] >. Region A of PKS 2349-014 correspond 
to a physical companion located 4" West of the nucleus and its 
emission-line ratios would place it in the composite region of 
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the BPT diagram. However, the H/3 and Ha lines are redshifted 
by ~200kms~' with respect to the [Om] lines. This can be ex- 
plained most likely due to an seeing induced superposition of 
emission lines corresponding to an H n region and AGN ionised 
gas and we classify this particular emission-line region as solely 
being due to star formation. These knots of star formation are 
likely to found in interacting/merging systems, besides a large 
AGN ionised region, as reported by Villar-Martm et al. ( 201 lb) > 
from observations of type 2 QSO at redshift 0.3-0.6. 

For the majority of regions we are only able to infer a single 
emission line ratio, either [Om] A5007/H/3 or [Nn] /16583/Ha, 
so that the ionisation mechanism can generally not be deter- 
mined unambiguously. We assumed that all regions with only an 
Ha detection are characteristic for H n regions. If these regions 
were instead ionised by AGN, the [O in] line would be at least as 
bright as the Ha line assuming a theoretical Balmer decrement 
of HafH/3 ~ 2.87 ( |Osterbrock & Ferland|2006) and a [Om]/H0 
line ratio of >3 as a threshold for AGN ionised gas. Dust ex- 
tinction along our line-of-sight within the host galaxy can dilute 
the [O m] and H/3 emission further, but a high dust extinction is 
often found around Hn regions and is directly proportional to 
the current star formation rate (SFR) (e.g. Hop kins et al.| 2001 ; 
|Sullivan et al.|20 01 ). Thus AGN ionisation would be an unlikely 
scenario in those cases. 

On the other hand, it is clear from the two branches in the 
BPT diagram that the [O m] A5QQlfH/3 line ratio alone is not suf- 
ficient to constrain the ionisation mechanism. An emission-line 
ratio of log([Om] A5001 /H/3) > 0.9 can still be securely as- 



signed to AGN ionisation which is the case for 7 objects. Ho 
|etaE1 ( [T997l > defined a threshold value of [Om] /15007/H/3 > 3 
above which AGN ionisation is likely to be the dominant pro- 
cess. This can be physically explained by the fact that the star 
forming branch in the BPT diagram represents actually a se- 
quence of increasing gas-phase metallicity from high to low 



[Oiii] /Hyg line ratios (e.g. |Dopita et aL]|2000| |Pettini & Pagel 
|2004] >. A high [Om] i5007/HyS line ratio can reasonably ex- 
cluded as a signature for star formation for luminous AGN 
by considering the mass-metallicity relation of galaxies (e.g. 
|Tremonti et al.||2004| l. Luminous AGN resid e almost exclu- 
sively in massive galaxies above M, > 1O 1O M (Heckman et al. 
|2004 , Kauffma nn et al.|2003| , corresponding to a high gas-phase 
metallicity (Z > Z Q ). Thus, we will rely on the Ho et al. criterion 
to classify the ionisation of our EELRs, but we caution that the 
argument as outlined above might not be correct for all EELRs 
as the gas-phase metallicity in galaxies typically decreases with 
distance from the galaxy centre (e.g. Henry & Worfhey|1999| l. 



From the detected regions of extended ionised gas we clas- 
sify 13 of 19 to be entirely ionised by the AGN, 3 of 19 to be 
powered by Hn regions only, and 3 of 19 exhibit signatures 
for both mechanisms dominating in spatially separated regions. 
Interestingly, all the H n regions are located in disc-dominated 
galaxies or in the strongly interacting system PKS 2349-014. 
|Letawe et al.| ( |2007[ l reported from their on-axis longslit spec- 
troscopy of luminous QSOs that the ISM was in most cases 
ionised either by young stars or shocks. They found a completely 
AGN-ionised ISM only in the case of bulge-dominated galaxies, 
while their disc-dominated system or morphological disturbed 
galaxies show signatures of Hn regions in almost every case. 
Since inactive disc-dominated galaxies are known to have young 
stellar population and continue to from star in their disc, it may 
not be surprising that this is true also for QSO host galaxies. 



5.4. Dependence on nuclear properties 

Compared to our earlier work presented in Hu08, we increased 
our sample by 1 1 objects and improved the QSO-host deblend- 
ing algorithm. Our main result in Hu08 was that the presence 
of an EELR is closely correlated with certain spectral proper- 
ties of the QSO itself, most prominently the EW of Fe n and the 
FWHM of the broad Hf3 line. Both are the main drivers for the 
Eigenvector 1 parameter space ( |Boroson & Green|1992| l describ- 
ing most of the variance among rest-frame optical QSO spectra. 
While this result is still valid, some of the new EELR detections 
require additional notes. 

We classified I Zw 1, RBS 219, PG 1700+518, 
PG 2130+099, and PHL 1811 as EELR non-detections in Hu08 
for which we now find evidence for an EELR. All of them have 
a narrower H/3 width than the proposed threshold at 4000 km s _1 
and significant Fen emission in their spectra. The main point 
about the EELRs around I Zw 1, RBS 219 and PHL 1811 is 
that all the extended emission is ionised solely by star forma- 
tion rather than AGN radiation. Thus, we refine our state- 
ment in Hu08 such that only the presence of extended AGN- 
ionised gas on kpc scales is correlated with the nuclear prop- 
erties as described above. Our improved deblending technique 
allows us also to detect a small ENLR close to the nucleus of 
PG 1700+518 and PG 2130+099. Similar to Mrk 1014, which 
we noted as an outlier in Hu08, these objects are known to host a 
triple radio source. Radio jets may therefore be an additional in- 
gredient altering the general trends, which we will explore more 
closely in the next sections. 

The additional 11 SDSS QSOs were specifically selected 
from their nuclear properties as candidates to show luminous 
ENLRs. This could be confirmed in most cases. Of the 4 EELR 
non-detections, 3 were observed through clouds (cf. Table |4]i 
which substantially reduced our chance to detect low surface 
brightness features. Thus, our statistics are not significantly in- 
creased compared to Hu08 and we focus on a detailed study of 
the ENLRP] The detected ENLRs in the remaining 7 new objects 
fit well into the correlations of Hu08, reinforcing the results re- 
ported in that paper. 



6. Extended narrow-line regions around QSOs 

6.1. The ENLR morphology 

As a first step, we classified the ENLR morphologically by eye 
into three categories: 'one-sided' (3 objects), 'biconical' (6 ob- 
jects) and 'round' (7 objects). Assuming elliptical isophotes, we 
also measured the ellipticity (1 - b/a) of the ENLR from the 
light distribution at a level of 2cr above the background noise 
in the reconstructed [Om] narrow-band images. Furthermore, 
we estimated the position angle (PA) of the major ENLR axis 
by searching for the maximum flux within bipolar cones of 20° 
opening angle. All these parameters are summarised in Table [8] 
The ellipticities of the ENLRs agree with our visual impression 
of round and biconical morphologies with a transition value at 
e ^ 0.4. An exception is PKS 2349-014 because its ENLR is 
quite lopsided and the southern high surface brightness region 
dominates the whole ENLR structure. 

We find a large fraction of 37.5% biconical ENLRs which 
even increases to 56% when the one-sided ENLRs are also as- 



4 Of course an ENLR is also an EELR, but we want to emphasize the 
AGN-dominated ionisation source with a different naming following the 
literature. Thus, an ENLR may either be embedded in a larger EELR 
structure, or representing the entire EELR. 
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Fig. 11. The nucleus-subtracted ENLR light distribution of the [O in] line in comparison to the resolved radio emission. For the RQQs the high 
spatial resolution VLA radio images are overplotted as blue contours, while the primary radio axis of the radio lobes in RLQs is indicated by the 
blue arrows. The white cross marks the position of the QSO on which the radio emission is centred. 



Table 8. Luminosity and morphological parameters of the ENLR. 



Name 


logi E NLR([Om]) 


Morph. 


e 


PA 




[ergs" 1 ] 






[°] 


detected ENLRs 


PG 0026+129 


41.80± 0.14 


round 


0.30 


-57 


PG 0052+251 


41.75 ±0.13 


one-sided 




34 


SDSS J0155-0857 


41.43 ±0.19 


round 


0.39 


5 


Mrk 1014 


42.22 ±0.14 


biconical 


0.60 


-65 


SDSS J0846+4426 


43.08 ± 0.09 


round 


0.09 


83 


SDSS Jl 131+2623 


42.11 ± 0.14 


biconical 


0.64 


18 


PKS 1217+023 


41.82 ±0.16 


biconical 


0.44 


62 


SDSS J1230+6621 


41.90 ±0.15 


round 


0.32 


5 


PKS 1545+210 


42.35 ±0.13 


biconical 


0.45 


-49 


PG 1612+261 


42.80 ± 0.12 


round 


0.20 


39 


SDSS J1655+2146 


43.03 ± 0.09 


round 


0.19 


43 


PG 1700+518 


41.81 ±0.18 


one-sided 




-27 


PG 2130+099 


40.90 ± 0.32 


one-sided 




-28 


HE 2158-0107 


42.32 ±0.13 


biconical 


0.46 


-46 


PKS 2349-014 


42.11 ±0.11 


biconical 


0.18 


6 


HE 2353-0420 


42.47 ± 0.12 


round 


0.15 


57 


undetected ENLRs 


IZw 1 


< 41.07 








SDSS J0057+1446 


< 41.90 








RBS 219 


< 41.57 








HE 0157-0406 


< 41.84 








SDSS J0948+4335 


< 42.66 








3C 273 


< 42.10 








HE 1228+0131 


< 41.51 








SDSS J1230+1 100 


< 42.86 








PG 1427+480 


< 42.42 








SDSS J1444+0633 


< 42.97 








HE 1453-0303 


< 41.77 








PHL 1811 


< 41.72 








RBS 1812 


< 41.51 








HE 2307-0254 


< 41.63 








PG 2214+139 


< 41.03 









sumed to be intrinsically biconical. The biconical fraction is in 
agreement with the one determined by Schmitt et al. (2003b) for 
type 1 Seyfert galaxies. However, it is expected that the open- 
ing angle of AGN ionisation cones increases with the luminos- 



ity of the AGN in the "receding torus" model ( |Lawrence|1991| l. 
Thus, biconical/elongated ENLRs should be less frequent than 
in lower luminosity Seyfert galaxies. This may indicate that the 
ENLR morphology is significantly affected by secondary prop- 
erties, such as the gas density distribution with respect to the 
ionisation cones (Mulchaey et al. 1996b| l, or by the presence and 
direction of a radio jet. 

A close alignment between the EELR and radio axis particu- 
larly around high-redshift RLQs and radio galaxies has been well 
established ( McCart hy" et al.|1995||1987| ), which i s often recog- 
nised also for lower luminosity Seyfert galaxies (Capetti et al.| 



1996 |Falcke et al. 1998;Naga~et al.|1999; Schmitt et al.p003b 



Wilson & Tsvetanov|1994| l on much smaller physical scales. In 
Fig. 11 we directly compare the ENLR light distribution against 
the radio morphology for those QSO with available interfero- 
metric radio data. Among the 7 objects for which the PA of the 
radio jet is known, only PKS 1545+210 is unaligned with re- 
spect to the radio axis, while the axes are aligned within ±20° 
(rms) for the majority of QSOs. This supports the view that ra- 
dio jets even with relatively low power (< 3 x 10 24 WHz" 1 at 
1 .4GHz) can have a significant influence on the morphology and 
characteristics of the ENLR. 



6.2. ENLR size measurements 

Previous ENLR surveys used various definitions for the ENLR 
size. Ground-based narrow-band images of 47 Seyfert galax- 
ies were analysed by Mulch aey et al.| ( |1996a[ ), who measured 
isophotal ENLR sizes at a limiting surface brightness o f 2 x 
10~ 16 ergs _1 cm -2 arcsec -2 . Circular annuli were used by |Ben-| 



Inert et al.| ( |2002| l to infer the ENLR radii at which the [O in] sur- 
face brightness dropped below the background surface bright- 
ness limit. Finally, Schmitt et al. (2003a) estimated the maxi- 
mum projected sizes of the [O m] light distribution by eye from 
their HST narrow-band images. Thus, different ENLR size defi- 
nitions are not directly comparable with each other, and some of 
them depend critically on the depth of the data. 

We measured three different ENLR sizes from the [O in] light 
distribution at 2<x above the background noise: 1) The rgs radius 
defined as the circular aperture containing 95% of the ENLR 
[Om] flux, 2) the effective radius r e , defined as the luminosity- 
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Table 9. Measured sizes of the ENLR. 



6.3. ENLR size-luminosity relation 



Name 


'"95 
["] 


i 

[' 


e 

'] 


' ISO 

["] 


rvj UUiUT li7 


3.4 = 


:0.7 


1 Q - 

1.7- 


C \J.J 


3.1 : 


:0.8 




3.4 = 


= 0.6 


1.7 - 


h 4 


4.1 : 


:0.8 


snss toi 55-0x57 


2.0: 


= 1.3 


1 - 
i .yj - 




2.0: 


: 1.5 


Mrk 1014 

1V11 Jv IVJlt 


> 3.9 




1 Q 


>4.3 


SDSS T0R4fS+4496 


2.6 = 


= 1.1 


1 ^ - 
i .j - 


h 4 


3.1 : 


: 1.4 


snss ti 1 3i +9693 


3.3: 


= 0.8 


1.4 - 


t 0.4 


4.5: 


: 1.0 


PICS 19174-0?^ 


3.3: 


= 0.6 


2.1 - 


h n 4 


3.2: 


:0.7 


snss H930+6691 


1.6 = 


= 0.9 


0.8 - 


i- 0.4 


2.0: 


= 0.8 


PKS 1545+210 


3.4: 


= 0.7 


1.8 = 


t0.4 


3.2: 


= 0.7 


PG 1612+261 


> 2.7 


> 


1.3 


> 3.3 


SDSS J1655+2146 


2.3: 


= 0.6 


1.0 = 


t0.3 


3.1 : 


= 0.7 


PG 1700+518 


3.6: 


= 1.0 


1.8 = 


t0.4 


2.2: 


= 1.8 


PG 2130+099 


2.5 : 


= 1.1 


1.3 = 


t0.5 


1.8: 


= 1.4 


HE 2158-0107 


> 3.4 


> 


1.9 


> 4.8 


PKS 2349-014 


3.1 : 


= 0.6 


1.8 = 


t0.3 


5.1 : 


= 1.2 


HE 2353-0420 


2.5: 


= 0.6 


1.1 = 


t0.3 


2.3: 


= 0.6 



weighted mean of projected pixel distances to the centre, and 3) 
the maximum isophotal radius r; so out to a limiting [Oiii] sur- 



' ISO 

face brightness of 2 x 1CT 16 erg s" 1 



cm 2 arcsec 2 . The measured 



ENLR sizes are reported in Table [9] with errors based on Monte 
Carlo simulations of the QSO-host deblending process. We gen- 
erated 200 different realisations of each object datacube consis- 
tent with the variance in each pixel, which we similarly analysed 
to infer lcr errors as the rms of the parameters. Because of the 
low spatial sampling of our reconstructed PMAS narrow-band 
images, we further added 0'.'5 to the uncertainty of rg$, and r lso 
as well as 0'.'25 to the uncertainty of r e . The seeing certainly in- 
creases the apparent size of the ENLR, but its strength strongly 
depends on the intrinsic emission line distribution and also on 
the size definition. We emphasise that any correction for the see- 
ing would be a small effect given the large apparent ENLR sizes 
of several arcseconds, so that we did not perform any seeing cor- 
rection. The ENLR apparently reaches beyond the PMAS FoV 
in three cases, HE 2158-0107, Mrk 1014, and PG 1612+261, 
for which we set lower limits to their ENLR sizes. For all other 
objects the surface brightness in the ENLR is sufficiently low at 
the edge of the PMAS FoV so that no significant extension of 
the ENLR is expected beyond the FoV. 

Because of the high frequency of elongated ENLR struc- 
tures, we tested by how much longslit spectroscopic observa- 
tions would underestimate the ENLR size if not properly aligned 
along the ENLR major axis. We extracted synthetic longslit sur- 
face brightness distributions from the narrow-band images along 
the major and minor axis of the ENLR. We find that ENLR sizes 
are underestimated on average by a factor of 2 when the slit is 
oriented along the major axis of the ENLR. For individual ob- 
jects even a factor of 6 can be reached. Thus, the choice of 
the longslit position angle has a stronger impact on the inferred 
ENLR size than the 30% reported by |Greene et al.|p0TT) based 
on their multiple longslit observations of type 2 QSOs. One 
possible explanation for this discrepancy could be the confusion 
of extended and unresolved [O m] emission of the nucleus, en- 
hanced by the seeing. We can deblend these two component in 
our QSO using the broad emission lines, but those are not visible 
in the spectra of type 2 QSOs by definition. 



In Fig. 12 we compare the projected physical ENLR size against 
the total [Oiii] luminosity (ENLR+QSO) (left panels) and the 
bolometric luminosity Lboi of the QSO (right panels). The ENLR 
size appears to be almost uncorrected with the [O in] luminos- 
ity given the low Spearman rank correlation coefficient and a 
high probability for no correlation (p = 0.27, P = 0.3, for r^). 
We find that the correlation coefficient is significantly higher for 
the QSO continuum luminosity L5100 0° = 0.65, P — 6 x 10" 3 , 
for rgs). The correlation between the ENLR size and the QSO 
continuum luminosity has not been investigated so far, because 
all previous studies focused either on type 2 AGN or on low- 
luminosity Seyfert nuclei which often exhibit a substantial con- 
tamination by the host galaxy continuum. |Bennert et al. (2002) 
compared the ENLR size also with the broad H/8 luminosity of 
the QSOs which is well correlated with the QSO continuum 
luminosity (e.g. |Greene & Ho||2005| |Punsly & Zliin^[20TT] 
Schu lze et al.|2009) >. Their data also imply a higher correlation 
coefficient for the ENLR size as a function of H/3 luminosity 
which however they did not notice in the paper. The fact that 
ENLR size is much better correlated with the AGN continuum 
luminosity than with the normally used [O in] luminosity shows 
that the [O in] luminosity is a less reliable proxy for the amount 
of ionising photons compared to the continuum luminosity, al- 
though the [Oiii] luminosity has been commonly found to scale 
well with the QSO luminosity (e.g. Heckman et al.|[2004 Za 
|kamska et al.|2003| >. 

Because reliable ENLR characterisation are scarce for lower 
luminosity Seyfert 1 galaxies, we cannot study the ENLR size- 
luminosity relation in detail in this paper and restrict ourselves 
to average quantities. The mean size of the detected ENLRs for 
our type 1 QSOs is r e = 5kpc ± 2kpc and = lOkpc ± 3kpc 
with a similar value for r; so at a median [Oiii] luminosity of 
log(L([Om])/[ergs -1 ]) = 42.7 ± 0.15. Just a few objects seem 
to follow the scaling relation of Bennert et al. ( 2002 ), while the 
majority are a factor of 2 larger than the extrapolation of their re- 
lation to the corresponding [O in] luminosities. The discrepancy 
is even larger for the scaling relation of Schmitt et al. ( 2003b| for 
their Seyfert 1 sub-sample. This is may be unsurprising given the 
fact that the limiting surface brightness of the HST narrow-band 



images is at least a factor of 10 higher. The samples of Greene 
|et al.| ( |201 1) and |Fu et aT7| ( |2012| > have [Om] luminosities and 
limiting surface brightnesses (~lxl0~ 16 erg s" 1 cirT 2 arcsec -2 ) 
comparable to our sample of type 1 QSOs. Their reported ENLR 
sizes are of the order of lOkpc, which matches with the results 
from our sample, although their ENLR sizes seem a bit larger at 



lower [Om] luminosities (cf. Fig. 12 1. However, when consider- 
ing the continuum luminosities for the few type 1 QSOs within 
the sample of |Fu et aT this discrepancy is largely reduced. 

The combined data of all studies strongly suggests that 
the ENLR scaling relation cannot be as different for type 1 
and type 2 AGN as implied by the results of Sch mitt et al.| 
(2003b). If we extrapolate their relation to an [Om] luminos- 
ity log(L( [Om])/ [ergs -1 ]) = 42.5, it would imply an ENLR 
size at least a factor of 2 lower for type 2 QSOs than for type 
1 QSOs. Additionally, it is obvious that the zero-point of the 
ENLR size-luminosity relation based on HST data is too low be- 
cause of the surface brightness sensitivity as discussed above. 
The slope of the relation is a different issue and may be a phys- 
ically more interesting parameter, but cannot be inferred from 
our small range in luminosities. Greene et al. ( |201 1[ ) determined 
a shallow ENLR size-luminosity relation only because ENLR 
sizes from Seyfert 2 galaxies were taken into account. Those 
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Fig. 12. ENLR sizes as a function of total [Om] luminosity (left panels) and the QSO continuum luminosity at 5100A (right panels). The used 
size definitions are isophotal radius r iso (upper panels), effective radius r c (middle panels), and the radius containing 95% of the ENLR flux r 95 
(lower panels). The solid, dashed and dotted lines in the lower left panel correspond to the ENLR size-luminosity relations of Bennert et al. ( 2002 ), 
Schmitt et al. (2003b), and Greene et al. (201 1 1, respectively. The black filled symbols represent RQQs whereas the red filled symbols highlight 
the few RLQs in our sample. The green opened circles are measurement for luminous type 2 QSOs as reported by Greene et al. (2011) and the 
blue triangles are measurements for double-peaked AGN (open symbols denote type 2 and filled type 1 AGN) provided by |Fu et al.H2012| . 



were estimated by extrapolation of the NLR surface brightness 



profiles to large radii based on longslit spectroscopy (Fraquelli 
|et al.|2003) . Whether this is a robust procedure to obtain reliable 
and unbiased (E)NLR sizes for low-luminosity AGN still needs 
to be confirmed by direct observations. 



|Netzer et al.| ( |2004[ l emphasised that the ENLR size- 
luminosity relation should not be a universal law. The ENLR 
size of luminous QSOs at high redshift, implied by the current 
relation, would exceed those of the host galaxies by more than 
an order of magnitude. Our QSOs are just at the limit where the 
ISM of the entire host galaxies is ionised, so that we do not yet 
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see a break in the ENLR size-luminosity relation. On the other 
hand, we were unable to detect any ENLR for about half of our 
QSO sample, at similar QSO continuum luminosities. Thus, it 
is not clear whether a universal ENLR size-luminosity relation 
does actually exist. A few of our ENLR non-detections are as- 
sociated with compact and apparently undisturbed host galaxies, 
such as PG 1427+480 (cf. Pig |7j, which might have a compara- 
bly small ENLR size that cannot be resolved by our observations. 
It is possible that the host galaxy size imposes a natural limit for 
the ENLR because of the sharp decline in gas density. Alterna- 
tively, galaxy interactions/mergers might increase the amount of 
available gas illuminated by the QSO in contrast to objects with 
non-detected ENLRs as suggested by Matsuoka] ( |2012[ ), which 
would enhance the visibility of the ENLR and possibly also its 
size. We already emphasized in Hu08 that the detection or non- 
detection of an ENLR correlates with certain spectral properties 
of the QSO that drive the variance along the Eigenvector 1 pa- 
rameter space (Bor oson & Green|1992 l. It is puzzling how the 
EW Fe ii in the nuclear spectrum can be causally connected with 
properties of the ENLR on kpc scales. We doubt that this can 
be ascribed solely by galaxy interactions; fundamental physical 
properties of the nucleus that are yet to be identified probably 
play an important role. 
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Fig. 13. [O in] luminosity of the ENLR as a function of the unresolved 
QSO [O m] luminosity. Squared symbols indicate our own PMAS mea- 
surements and star symbols refer to the literature data of Stoc kton "&] 
MacKenty (1987). Radio-quiet and radio-loud objects are marked as 
black and red symbols, respectively. The open grey symbol represent 
upper limits and the dashed line is a linear fit to our PMAS data (squared 
symbols) only. 



6.4. ENLR luminosities 

In addition to the [O m] luminosities of the resolved ENLR we 
also obtained robust upper limits for Lenlr([Oiii]) in the unde- 
tected cases. We added ENLRs with the observed shapes to the 
original datacubes until the input luminosity reached the detec- 
tion limit. We consider an ENLR to be detected at 3<x signifi- 
cance when the fluxes in 9 adjacent spaxel were above the 2<x 
background noise level in 95% of 200 different Monte Carlo re- 
alisations. Here we implicitly assume that the ENLR luminosity 
is the limiting factor to detect an ENLR and not its size. The 
[Om] luminosity of the ENLRs (Lenlr([Oiii])) is reported for 
all QSOs in Table [8] and compared with the unresolved [O in] 



luminosity of the nucleus (Lqso([Oiii])) in Fig. 13 Our mea- 
surements indicate a correlation of the two different luminosi- 
ties that is in agreement with the derived upper limits. With the 
Astronomy Survival Analysis package (AS URV, |La valley et al.| 
1992| l we determine a generalised Spearman rank correlation co- 
efficients of p = 0.76, taking the censored data into account. The 
probability for a non-correlation is P < 10~ 5 for the combined 
data set, which covers more than 2 orders of magnitude in [O in] 
luminosity. This possibly suggests that the nebulae are the ex- 
tensions of the unresolved NLRs in terms of the covering factor 
for the ionising radiation of the AGN. A linear regression analy- 
sis with the Buckley-James method as part of ASURV yields a 
relation of 

log(L ENL R([Om])) = (1.7±0.3)log(L QSO ([Om]))-(30.4±2.2) 

(5) 

A similar trend was already reported from the [O ni] narrow- 
band survey of low-redshift QSOs by |Stockton & MacKenty| 
( |1987| l and the narrow-band survey of local Seyfert galaxies car- 
ried out by Mulchaey et al. ( 1996a). The detected EELRs in the 
|Stockton~& MacKenty ( 1987 ) sample are mostly associated with 
RLQs and appear to follow a similar scaling relation as implied 
by our data at high [Om] luminosities. That their measurements 
are systematically lower than our relation can be easily explained 
by the fact that they excluded a central aperture with a radius of 
~3" to avoid contamination by the QSO light. 
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Fig. 14. [O m] luminosity of the ENLR as a function of QSO contin- 
uum luminosity at 5100A. Symbols are as defined in Fig. 13 



If the [Om] luminosity was in general an indicator for the 
amount of ionising photons emitted by the AGN, we would also 
expect a significant correlation of the ENLR luminosity with 
the QSO continuum luminosity L 51 oo- We already reported in 
Hu08 that this is not the case and present the corresponding dia- 
gram now for our extended sample in Fig. [14] The generalised 
Spearman rank correlation coefficient computed with ASURV 
is p = 0.108 with a probability of P = 0.55 to be completely 
uncorrected. This is rather surprising because of the strong cor- 
relation between the ENLR size and the continuum luminosity 
as shown in Fig. 



12 



To investigate the cause of this lack of correlation, Mat suoka] 
(2012) compared the AGN continuum luminosity against the Ed- 
dington ratio using a large compilation of available QSO data in- 
cluding the one presented in this article. They quote to see a clear 
anti-correlation between the Eddington ratio and the ENLR de- 
tection rate and argue that this is likely caused by the availability 
of gas given the signs of interaction found in these systems (e.g. 
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Fig. 15. Distribution of QSOs with detected and non-detected ENLRs 
in the M Bll -L bo i plane. The ENLR luminosity or its upper limit is en- 
coded in the symbol size. QSO with detected ENLR are indicated by 
filled symbols and non-detections by grey open symbols, respectively. 
Solid and dashed lines correspond to constant Eddington ratios of A = 1 
and A = 0. 1 , respectively. 



|Husemann etlQT1|20"T0l |Villar-Martm et d.pOTTb) . We show 
the distribution of our QSOs with ENLR detections and non- 
detection in the Mbh-L\, \ plane in Fig. 15 Many QSOs with a 



non-detected ENLR appear to populate a similar region than the 
ones with a detected ENLR. A two dimensional Kolmogorov- 
Smirnov (K-S) test (Fasano & Franceschini 1987 ; Peacock 1983) 
confirms a high probability of P = 0.74 that both samples 
(ENLR detections and non-detections) are drawn from the same 
parent population, which questions the trend proposed by |Mat-| 
suoka]( [2012| l. What we can confirm is that we see a deficit of 
ENLR for QSOs close to Eddington ratio (log A > -0.5), but 
the number of QSOs in our sample is too small in that regime to 
make a robust statements. 



Nevertheless, the ENLR luminosity apparently cannot be de- 
termined solely by the amount of ionising photons, so that it 
seems indeed required to invoke the mass and distribution of gas 
in the host galaxy as additional parameters. Another aspect is the 



increasing evidence supporting the clumpy torus model ( Alonso 



|Herrero et al.|201 l[|MoFet al.|2009||N"enkova et al.|2008| >. In this 
model it is just a matter of probability if our line-of-sight is inter- 
cept by dusty gas clouds of sufficient column density to obscure 
the nucleus. Thus, the intrinsic fraction of radiation that is able 
to escape into the host galaxy can have a large range even for 
apparent type 1 QSOs. It is therefore not surprising that the cov- 
ering factor has been shown to be the dominant parameter that 
controls the [Om] luminosity in QSO spectra (Baskin & Laor| 
|2005] >. Interestingly, it was recently reported by |Xu et aL ( 2012[ l 
that also the gas density of the ISM correlates with the nuclear 
properties in the sense that objects with narrow H/3 lines and 
stronger Fe n emission have a systematically denser ISM. This 
could imply that our non-detected ENLRs may be explained sim- 
ply by the low efficiency of ionising photons to ionise the dense 
ISM clumps due to self-shielding, It may also even explain the 
presence of star formation, as high gas density clumps are more 
favourable to have cool molecular gas cores. 
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Fig. 16. [O m] luminosity of the ENLR against radio luminosity at 
1.4 GHz. Black solid symbols refer to objects with detected and grey 
open symbols to objects with undetected ENLR (upper limits). 



6.5. Dependence on radio properties 

Radio observations have shown that the nuclear [O in] luminos- 
ity is correlated with the continuum radio luminosity in radio 
galaxies and RLQs (e.g. Baum & Heckman 1989] [Rawlings & 



Saunders 



Ulvestad| 



l99T)|Xu et al.|1999|> as well as in RQQs (e.g. |Ho &" 
2001] [Miller et al.||19931 |Xu et al.||1999> . The radio 



(6 cm) to optical (4400 A) flux ratio, the R parameter, has of- 
ten been used to distinguish between RQQs (R < 1) and RLQs 
{R > 10), so that their reported [Om] scaling relations often 
appear significantly offset from each other in terms of radio lu- 
minosity. Because the radio emission in RLQs is often dom- 
inated by their extended lobes, the two classes may follow a 
common relation only if the radio core power is considered ( |Xu| 
et al.||1999] l. Here we compare Lenlr([Oiii]) with the total ra- 
dio luminosity at 1.4 GHz (Fig. [To*] excluding the few known 
RLQs in our sample. It is difficult to argue about a possible cor- 
relation between these two quantities because of the numerous 
upper limits in both axes, radio and ENLR luminosity. The gen- 
eralised Spearman rank correlation coefficient is p = 0.56 with 
P = 0.003 for the non-correlation case. Nevertheless, it seems 
that the ENLR luminosity depends more strongly on the radio 
luminosity than on the QSO continuum luminosity, which we 
find is supported by three facts: 

1. An ENLR is detected around all QSOs with 
log^l^GHzj/tWHz" 1 ]) > 23.6 (8 objects). 

2. Ten of thirteen RQQs with a detected ENLR (77%) have a 
radio detection, while only five of fourteen RQQs with a non- 
detected ENLR (40%) are detected in the radio. 

3. In terms of radio-loudness, we detected an ENLR only 
around 26% of the RQQs with R < 1, but around 87.5% 
of the QSOs with 1 < R < 10. 

The ENLR luminosities of RLQs (not included in Fig. [16} 
are of the order of their radio-quiet counterparts, although the 
integrated radio luminosity is more than two orders of magnitude 
higher. Whether this is an intrinsic property or caused by the 
limited FoV of our IFU observation is something we are not able 
to say at present. 

|Sulentic et al.| ( |2000] > found that the broad H/3 line 
width of RLQs is almost always >4000kms _1 , and de- 
fined two different populations of QSOs. Population A with 
FWHMh^ <4000kms _1 is entirely dominated by RQQs and 
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log<ai.4GHz)/[WHz ]) 

Fig. 17. ENLR size as a function of radio luminosity at 1.4 GHz. 
Measurements for our RQQs with detected ENLR are indicated by the 
black filled symbols and open green symbols denote measurements for 
type 2 QSOs as reported by Gre ene et al.| ( [201 1) . 



Population B with FWHM^ >4000kms _1 contains all RLQs 
and a number of RQQs with otherwise similar characteristics. 
These two populations represent opposite regimes along the 
famous Eigenvector 1 (Bor oson & Green] |1992[ ), which ulti- 
mately links the width of H/?, the strength of the Fen com- 
plexes and the relative strength of [Om] with respect to H/3. 
The physical drivers for this sequence in parameter space are 



thought to be Mbh and A (Boroson 
2003) 



2002 Marziani et al. 



uminous ENLR is pref- 



In Hu08 we emphasised that a 
erentially found around QSOs with FWHM^g > 4000 km s~' 
and weak Fen emission that would correspond exclusively to 
Population A. Our extended sample now contains a few ob- 
jects, such as the radio-intermediate QSOs PG 1612+261 and 
SDSS 1655+2146, which are surrounded by a luminous ENLR 
but have a FWHM^ < 3000 km s _I and therefore do not cor- 
respond to Population A as one would have naively expected. 
This is in agreement with the results of Za mfir et al.| (|2008 ) who 
found that radio-intermediate QSOs also cover the full range in 
broad H/3 widths like RQQs and are not restricted to Population 
A like RLQs. 

Furthermore, Leips ki et al.| ( |2006[ l reported a significant 
trend between the size of the radio jet in RQQs and the ENLR. 
Because the radio emission is resolved only for very few of our 
QSOs, we investigate instead whether the ENLR size depends 
on the radio luminosity (Fig. 17 1 including the recent measure- 
ments for type 2 QSOs by |Greene et al.| ( |201 l| l. Only a weak 
trend between the two quantities seems to be present if there is 
one at all. Given our tentative evidence that the radio luminos- 
ity is more strongly correlated with the ENLR luminosity than 
its size, we speculate that a radio jet increases the covering fac- 
tor of the gas for the ionising photons of the nucleus, reaching 
out to larger distances. This would well match with the scenario 
proposed by Hopki ns~et al.| ( |20T0) >, who argued that a weak out- 
flow in a hot diffuse medium can possibly increase the impact 
of the AGN on the host galaxy via a two-stage process. First, 
cold molecular gas clouds expand or are destroyed by the out- 
flow induced pressure gradients and instabilities, which could be 
generated by a radio jet. Second, the molecular clouds are more 
efficiently illuminated by the QSO radiation, because of their in- 
creased cross-sections and covering factor, so that the AGN can 
dissipate more energy into the medium. 



7. Kinematics of the extended ionised gas 

7.1. Construction of 2D velocity fields 

In order to study the kinematics of the EELR we used the bright- 
est emission lines to infer the 2D radial velocity distribution of 
the ionised gas. We modelled the [Om] A/14960, 5007 doublet 
and/or the Ha emission line adopting a simple Gaussian line pro- 
file in each spaxel of the nucleus-subtracted datacubes. The line 
widths of the doublet lines were coupled and their line ratio was 
fixed to the theoretical value. We also modelled each spectrum 
assuming a Gauss-Hermite line profile up to the fourth order fol- 
lowing |van der Marel & Franx| ( |1993[ ) to identify asymmetric 
line profiles in the EELR not caused by noise. Asymmetric lines 
are not frequently found in the EELR of our QSO sample (see 
co-added spectra in Fig. [7j. In the following, we first present 
the 2D velocity fields for the bulk gas motion only. When line 
asymmetries are significant we replace the radial velocity of the 
single Gaussian fit with the mean radial velocity of the Gauss- 
Hermite result at that location. In addition to the 2D velocity 
maps, we constructed synthetic "longslit" curves along a certain 
PA directly from the maps in a geometrically manner. We aver- 
aged the velocities of all spaxel that are covered by a hypothet- 
ical slit of 1'.'3 width at an angular separation of 0'.'5 along the 
slit axis. The QSOs PG 1700+518, PHL 181 1 and PG 2130-099 
were excluded from this 2D analysis, as their EELR are too small 
in spatial extent. Still, it is worthwhile to note that we measured 
a redshifted radial velocity of ~800kms~' for the EELR struc- 
ture of PG 1700+518 (cf. Table|7|. Furthermore, we needed to 
exclude RBS 219 and SDSS J0155-0857, for which the low S/N 
of the extended emission lines did not allow a spatially resolved 
analysis. 

We used the narrow [O m] line in the QSO spectrum to es- 
timate the systemic redshifts of the QSO host galaxies. Previ- 
ous studies showed that the spatially unresolved [Om] line is 
often systematically blueshifted in QSOs (e.g. 



Boroson 



Eracleous & Halpern|2004l |Komossa & Xu|2007| |ViTK-M artfn 



2005 



et al.||20lTa I, which has been interpreted as a signature for cir- 
cumnuclear outflows. Whether such outflows commonly ex- 
tend out to kpc scales is one of the issues we want to address. 
We present the 2D velocity fields in three categories, (i) appar- 



ently large scale rotational motions (Fig. 18 i, (ii) velocity fields 



of morphologically identified major mergers (Fig. [T9j), and (iii) 
miscellaneous velocity fields (Fig. [2D) , which we discuss sepa- 
rately in the following. 



7.1.1. EELRs dominated by rotational motion 

We find signatures of large scale rotational motion around 
I Zw 1, PG 0052+251, PKS 1545+21, PG 1612+261, 
HE 2158-0107 and HE 2353-0420. The corresponding gas ve- 
locity fields are shown in Fig. [18] Because the turn-over point in 
the velocity curve is not reached in most cases, we modelled the 
synthetic longslit curves with a straight line (v r = Av r ■ r + vq) 
to search for the major kinematic axis with the steepest velocity 
gradient Av r among all -90° < PA < 90° in steps of 3°. The 
velocity curves along this kinematic major axis are also shown 
in Fig. [18] and the corresponding parameters are given in Ta- 
The majority of these systems (4/6) are morphological 



ble 10 



disc-dominated galaxies for which the kinematic major axis is 
well aligned with the photometric major axis within ±14° (rms). 
Thus, we suggest that the ionised gas is following the gravita- 
tional motion of the stellar disc in these cases. 
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Fig. 18. 2D velocity maps for objects with apparently large-scale rotational motion gas motion. The continuum light distribution of broad- 
band images are overplotted as thin solid black contours for comparison. A "longslit" velocity curves along the apparent major kinematic axis 
is extracted from the maps as described in the text. A thick solid black line indicates the orientation of the major kinematic axis. The radial 
velocity are estimated for each spaxel by modelling the corresponding emission line with a Gaussian profile. The vertical dashed line highlights 
the position of the QSO along the synthetic slit and the horizontal dashed line indicates the measured velocity zero-point of the velocity curve. 



To verify the reliability of our velocity fields, considering 
that the bright QSO has been subtracted before, we compare our 
ionised gas kinematics of I Zw 1 with 12 CO(1-0) data on scales 
of 1"— 14" ( |Schinnerer et al.|19 98 ). A radial Ha velocity of v r - 
166 + 22kms~ 1 at a distance of 2" is in goo d agreement with 
a velocity of 180 kms" 1 at 1.6" reported by |Schinnerer et al. 



( 1998 ) from the CO kinematics. Such a comparison between CO 
and ionised gas kinematics is only possible for I Zw 1, because 
CO observations with the required spatial resolution and depth 
do not yet exist for other QSOs in our sample. 



PKS 1545+210 and possibly also HE 2158-0107 are bulge- 
dominated host galaxies for which the inclinations of the puta- 
tive rotating gas discs are unknown. It is unlikely that the gas 
kinematics are caused by AGN-driven bipolar outflows, because 
the radio jet of PKS 1545+210 is oriented perpendicular to the 
kinematic major axis (see section [7T2| >, and HE 2158-0107 is 
a RQQ with a low radio luminosity. Interestingly, an apparent 
companion galaxy is located ~3" to the West of PKS 1545+210 
in the direction of the kinematic major axis and HE 2158-0107 
features an extended tail of gas . Both case are suggestive for the 
accretion of gas from companions or the environment itself (see 
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Table 10. Kinematic properties of rotational velocity fields. 



Name 


PA" 


Av r b 


V0 C 


R d 


Vr(Rf 




[°] 


[^] 




["Ape] 


[*P] 


IZw 1 


-28 


88 


704 


3/3.5 


166 


PG 0052+251 


2 


94 


134 


3/8.1 


284 


PKS 1545+210 


-46 


165 


3 


2/8.2 


299 


PG 1612+261 


6 


57 


39 


3/7.0 


119 


HE 2158-0107 


-49 


17 


5 


2/6.9 


49 


HE 2353-0420 


20 


70 


-25 


2/7.3 


138 



Notes. (a) Position angle along the highest velocity gradient. <b) Radial 
velocity gradient along the major kinematic axis per arcseconds. (c) Ra- 
dial velocity zero point that the QSO position with respect to the rest- 
frame defined by the nuclear [O in] line. <d) Characteristic radius from 
the centre used for a dynamical mass estimate <e) Radial velocity at the 
radius given given in (d). 



the case of HE 2158-0107 in |Husemann et al.|201 1| >, which will 
govern the rotation axis of the gas. 

A direct comparison of the Ho- and [O in] kinematics is only 
possible for PG 0052+251. We find that the [Om] kinematics 
does not match the clear rotational Ha velocity field with an off- 
set of up to ~100kms _1 . This kinematic separation supports our 
interpretation of the emission-line ratios in section 5.3 in which 
the Ha emission originates dominantly from H n regions in the 
discs, whereas the [Om] emission more often represents AGN 
ionisation cones. 



Our PMAS observation of Mrk 1014 covers only the inner 
lOkpc part of the QSO host (cf. Fig. [7}, so that we cannot trace 
the bulk motion of the ionised gas across the whole galaxy. The 
[Om] velocity curve along the EELR major axis appears to be 
peculiar. The symmetric curve has two sharp peaks in the radial 
velocity at +2" distance from the nucleus at which the radial 
motion changes from receding to approaching and vice versa. 
These features are associated with a region of strong jet-cloud 
interactions as discussed below in section [772] likely leading to 
non-relaxed bulk motion of the gas. 

The EELR of SDSS J 113 1+2632 is very elongated, and 
we can only study the kinematics in a part of the galaxy sim- 
ilar to Mrk 1014. The most striking feature is a giant shell 
or arc of ionised gas, 4" South-East of the nucleus, offset by 
352 + 33kms _1 from the systemic velocity in the rest-frame. 
This is a unique feature among all the QSO in our sample and 
does not correspond to any continuum substructure in the broad- 
band image of the complex merging system. One possibility is 
that the gas belongs to a faint companion galaxy illuminated by 
the AGN radiation. 

A strong velocity gradient is detected in the ionised gas of 
PKS 2349-014. The main axis of this gradient follows the large 
curved wisps that were identified in the HST images by |Bah-| 
call et al. ( 1995]>. Thus, the gas motion is clearly driven by 



gravitation. However, the radio-jet axis of this RLQ is oriented 
North-South, and the radio jet seems to kinematically disturb the 
ionised gas at the interception point of jet and ENLR. 



7.1 .2. Velocity fields of interacting QSO host galaxies 

In Fig. [19] we present the gas velocity fields of Mrk 1014, SDSS 
Jl 131+2632 and PKS 2349-012, which we identified to be on- 
going major mergers based on their optical continuum images. 
Because of the complex mass distribution and tidal forces, we do 
not necessarily expect a clear rotational velocity pattern in these 
systems. 



7.1.3. Miscellaneous velocity fields 

The radial velocity fields of the remaining four QSO hosts, 
PG 0026+129, SDSS J0836+4426, PKS 1217+023 and 
SDSS J1655+2146, are either asymmetric or lack a significant 
velocity gradient across the galaxy. All four host galaxies ap- 
pear to be bulge-dominated, in which the gas could be kinemat- 
ically hot and slowly rotating. However, the SDSS images of 
J1655+2146 and J0836+4426 do not rule out a disc component, 
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which would need to be studied with higher spatial resolution 
imaging. We do not find any strong evidence for extreme kine- 
matic features that would indicate non-gravitational motion. 



7.2. Evidence for jet-cloud interactions on kpc scales 

We already discussed a close connection between the ENLR and 
the presence and morphology of radio jets emerging from the 
nucleus in Section [6] To follow up on this we now address the 
question whether the kinematics of the ionised gas is also signif- 
icantly influenced by the jets. 

We detect significantly asymmetric [O in] line profiles close 
to the nucleus of Mrk 1014 and PG 1612+261 right at the loca- 
tion of their known radio hot spots (Fig. 21 1. The line width and 
radial velocity offset are 508 ± 84 km s _1 and 
for Mrk 1014 and 424±24km s" ' 



1012 ±48 km s" 1 
1 (FWHM) and -220±20kms" 1 
for PG 1612+261 at region A, respectively, indicating clearly 
non-gravitational motions. In addition to the blue-shifted com- 
ponent in Mrk 1014, a broad redshifted component is blended 
with the [Om] narrow component. The extended ionised gas 
kinematics of Mrk 1014 was already studied in detail by means 
of long-slit fLeipski & Bennert||2006[ ) and IFU spectroscopy 
( |Fu & Stockton |2009j Our data for this object are fully 
supporting their results. In addition, SDSS J0836+4426 and 
SDSS J1655+2146 similarly show asymmetric [Om] line pro- 
files close to the nucleus, which we separated into two differ- 
ent components. It is possible that jet-cloud interactions play 
a significant role in producing these complex kinematic fea- 
ture as well given their significant integrated radio luminosity of 
2 x 1 24 W Hz 1 and 4 x 1 23 W Hz 1 , respectively. Yet, no inter- 
ferometric radio images are available for these two objects that 
would be needed to confirm a clear association between radio 
emission and line asymmetries for a jet-cloud interaction sce- 
nario. 

The EELR of PG 1700+518 is almost coincident with the 
south-east radio hot spot, considering the low spatial resolution 
of the optical compared to the radio observation. It is tempting to 



speculate that the radio-jet is accelerating the gas to the observed 
velocities (~ 800 kms -1 ). On the other hand, the EELR may also 
belong to a companion galaxy or a starburst region as part of the 
merger, because a knot in the H-band adaptive optics image was 
identified by Marquez & Petitjean (2003 ). The brightest region 
of the EELR around the RQQ PG 2130+099 is matching with 
the southern hot spot of the triple radio source. However, we 
could not find any kinematic impact on the gas as in the other 
cases. 

The [Om] kinematics of the major merger host of the RLQ 
PKS 2349-012 are mainly quiescent. However, we find red- 
shifted [Om] emission by 400 ± 70km s with an intrinsic line 
width of ~700 km s (FWHM) at the interception point with the 
radio jet axis (2" North of the nucleus). This highlights the im- 
pact of the radio jet on the kinematics at the interception point, 
but not in the surrounding regions. 

The RLQs PKS 1217+023 and PKS 1545+210 reside both 
in bulge-dominated galaxies. The jet axis of PKS 1217+023 is 
aligned with the EELR major axis along which the velocity curve 
does clearly not agree with rotational pattern. Whether the gas 
is directly interacting with the jet remains uncertain because of 
the narrow line width and radial velocity <200kms _I . An inter- 
action of the jet with the EELR around PKS 1545+210 is most 
unlikely because the radio jet is aligned with the rotation axis of 
the gas disc. Such rotating gas discs are common in low-redshift 
radio galaxies (e.g. |Baum et al.||1992{ |Kotanyi & Ekers||1979| l 
with Centaurus A being the most prominent example. 

To summarize these individual results, in Fig. [22] we com- 
pare the maximum radial [O m] EELR velocity offset with re- 
spect to the zero-point of the velocity curves against the radio 
luminosity for all RQQs in our sample and that of |Greene et al.| 
(2011). A tentative trend is emerging that the radial velocity of 
the highly ionised gas on kpc scales increases with the radio lu- 
minosity, i.e. jet power, even in RQQs. Despite some contamina- 
tion with purely rotational motions at low radio luminosity and 
the unknown geometry and inclination of the gas motion, this 
supports the notion of a the systematic influence of radio jets on 
the ionised gas kinematics on host galaxy scales. A quite sim- 
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Fig. 21. Left panel: Strength of the third Gauss-Hermite coefficient h 3 representing an asymmetry in the line profile found to be significant 
only in 4 galaxies, i.e. Mrk 1014, PG 1612+261, SDSS J1655+2146, and SDSS J0836+4426. The black contours represent the [Om] surface 
brightness distribution for comparison and the blue cross marks the QSO position for reference. Right panel: Modelling of [O m] /14960, 5007 
doublet line in the co-added spectra (of the spaxels in the magenta box on the left) with two Gaussian components. The red line is the model to 
the measured spectrum shown in black, and the green line represents the individual Gaussian components. 



ilar trend was already reported by Miiller-Sanc hez et al.| ( |20ll) 
for sub-kpc scales in lower luminosity Seyfert galaxies based 
on near-infrared IFS. Nevertheless, our inferred radial velocities 
are in the majority of cases below th e escape velocity, w hich is 
of the order of 500-1000 km s" 1 (see |Greene et al.|201l) . Such 
a trend with increasing radio luminosity may be expected as the 
jets will carry mechanical energy and momentum of which a cer- 
tain fraction will certainly be dissipated in the surrounding gas. 
Radio-loud AGN are therefore known to exhibit jet-cloud inter- 



actions that accelerate the gas to significant velocities (e.g.LBald 
win et al.|1980l|Christensen et al . 12006 UHolt et al ||20031 |Mor- 



ganti et ai.||2005| |Nesvadba et al.||2006| [Sanchez et al.| |2004a>. 



On the other hand, detailed studies of some RQQs (e.g. SDSS 
J1356+1026,|Greene et al.||20T2] > and ULIRGs (e.g. Mrk 231, 
Rupke & Veill eux|201 1) have shown that a radio jet alone is en- 
ergetically insufficient to power high-velocity outflows present 



in these objects. The lack of high velocity ionised gas clouds 
in the majority of luminous RQQs is not necessarily in contrast 
with those observations. It could be that radio jets dominate the 
cloud acceleration even energetically in the majority of cases, 
but with some particular exceptions. An alternative explanation 
could be that radio jets significantly increase the effective energy 
dissipation of the QSO radiation into the ambient medium. We 
mentioned already in Sect. |6.5| the possibility that a jet might de- 
crease the dumpiness of the surrounding medium so the QSO 
radiation can act more strongly as proposed by Hopkins et al. 
(2010) in their two-stage AGN feedback model. 
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Fig. 22. Extended maximum [OIII] radial velocity in the ENLR of 
RQQs as a function of radio luminosity at 1.4 GHz. Only objects with 
detected radio emission are included in the diagram. Data from our IFU 
study of RQQs are indicated by the black filled symbols and open green 
symbols are from the long-slit data of type 2 QSOs (Greene st al.|201l) . 



7.3. The incidence of outflows among RQQs: No evidence 
for AGN feedback? 

AGN feedback has generally become an important ingredient in 
cosmological simulations of galaxy formation to efficiently stall 
star formation on short time scales. In these models either a sig- 
nificant amount of gas is expelled from the host galaxy by AGN 



outflows (e.g. Hopkins et al. 2005 ; Sanders et al. 1988) or the gas 
in the galaxy halo is episodically heated by radio jets to prevent 
further condensation of cold gas (e.g. Best et al. 2006 Bower 
|et al.|2006{|Croton et al.|2006) >. Clear observational evidence for 
the latter process are the intra-cluster X-ray cavities driven by ra- 
dio jets of the brightest cluster galaxies (see |Cattaneo et al.] 2009 ; 
Fabian 2012, and references therein). On the other hand, obser- 
vational signatures for AGN driven outflows are still scarce, and 
their interpretations are a subject of debate. 

Galactic superwinds with outflow velocities of up to 
1000 km s~' in the ionised gas have been detected in Ultra Lu- 



minous Infrared Galaxies (ULIRGs) (e.g. Harrison et al 



2012 



Heckman et al.|2000|IRupke et al.|2005j|Rupke & Veilleux|201 1 
Westmoquette et aL[ 2012i and massive post-starburst galax 
ies (Tremonti et al. 2007). These extraordinary outflows may 



be driven by the radiation of the central AGN with a com- 
parably low contribution from the starburst. Whether AGN- 
driven outflows are solely responsible for such outflow velocities 
is however questioned by recent observation of post-starburst 
galaxies that display high-velocity outflows without clear evi- 
dence for AGN (Diamo nd-Stanic et al.||2012| ). This is consis- 
tent with a comparison study of the radial gas velocity in X-ray 



selected AGN hosts and post-starburst galaxies by Coil et al. 
( |201 l| l. They concluded that while typical outflow velocities 
of 200 km s -1 are seen in both samples, the winds are likely 
driven by vigorous star formation and not by AGN. Furthermore, 
ionised gas with radial velocities exceeding more than 500km s _I 
have been detected mainly around radio-loud AGN (e.g. |Fu &| 
|Stockton||2009l|Holt et al.||2008] |Nesvadba et al.|[2006| , which 
can be attributed to the mechanical energy of jet. 

In our sample of 3 1 low-redshift QSOs we find signatures of 
outflowing (> 400 km s _1 ) ionised gas on kpc scales only in 3 ob- 
jects. In all these cases we identify a radio jet that is most likely 
driving the outflow. Quiescent gas kinematics, or no ionised gas 



ionised gas of all but one type 2 QSO, but they attempt to use the 
wings of the [O in] line to argue for the presence of some level 
of AGN feedback. The low outflow detection rates agree with 
the result of Krug et al. (2010), who compared the radial veloc- 
ity of the Na i D absorption line in the spectra of infrared-faint 
(SFR<10M G yr _I ) Seyfert galaxies at z < 0.05 and reported an 
outflow detection rate of 6% for Seyfert 1 and 18% for Seyfert 2 
galaxies, compared to more than 50% for starburst galaxies. 

Our observations clearly show a lack of ionised gas outflows 
in radio-quiet AGN at low redshift, which may challenge the in- 
corporation of AGN feedback in some theoretical models. How- 
ever, no ionised gas at all is found around almost half of the 
QSOs in our sample. These QSOs typically show weak [Om] 
emission in their nuclear spectra, which could be interpreted in 
two ways: (i) The majority of gas has already been expelled from 
the host galaxy (e.g. |Netzer et al.|2004| l, (ii) the majority of gas 
within the host galaxy is not sufficiently ionised by the AGN, or 
(iii) the galaxies are intrinsically poor of gas. With the advent of 
new radio/sub-mm telescopes, like the Atacama Large Millime- 
tre Array (ALMA), we have the opportunity to study the amount, 
distribution, and kinematics of the cold molecular gas at a suf- 
ficient spatial and spectral resolution. Given that the molecular 
gas is the seed of star formation, it may be the more direct tracer 
for AGN feedback processes than the ionised gas phase. 

8. Conclusions 

In this paper we presented a detailed analysis of the extended 
ionised gas around a large sample of 3 1 low-redshift QSOs ob- 
served with the PMAS integral field spectrograph on the 3.5 m 
telescope at the Calar Alto observatory. This is currently the 
largest sample of low-redshift type 1 QSOs observed with in- 
tegral field spectroscopy in the optical wavelength regime. In 
Husemann et al. (2008) we already studied the presence of ex- 
tended ionised gas in relation to the spectral properties of the 
QSO nucleus based on a subsample already observed at that 
time. With the full sample now in hand, we investigated the 
ionisation states, luminosities, sizes of the extended narrow line 
regions (ENLRs), and the kinematics of the ionised gas. In or- 
der to separate the unresolved QSO emission from extended host 
galaxy emission in the IFU data we introduced our new software 
tool QDeblend 3D . The main results of our analysis can be sum- 
marised as follows: 

- We detect extended ionised gas around 19 of 31 QSOs 
(61%), of which 13 are ionised mainly by the QSO (the 
ENLR), 3 are Hn-like regions produced by ongoing star 
formation, and the remaining 3 objects show signatures of 
both types in different extended regions. 

- Previous ENLR sizes obtained from HST narrow-band 
images are underestimated by at least a factor of 2 because 
of the much shallower surface brightness limits. Nine of 
sixteen ENLRs in our sample are elongated which means 
that longslit observations will strongly depend on the slit 
orientation. In the worst case they will underestimate the 
size by up to a factor of 6 for individual object, and a factor 
of 2 in the mean. 

- The typical size of the ENLR is rgs = 10 + 3 kpc at a median 
[Om] luminosity log^tOmD/tergs^ 1 ]) = 42.7 + 0.15 
of our QSO sample. Recent studies of type 2 QSOs at 



Article number, page 28 of 



34 



B. Husemann et al.: The properties of the extended warm ionised gas around low-redshift QSOs 



comparable luminosity by |Greene et al.| ( [201 l| l and |Fu et al.| 
" 201 2} obtained very similar ENLR sizes compared to our 



type 1 QSOs. This implies that the ENLR size-luminosity 
relation for type 1 and type 2 AGN cannot be as different 
as implied by the results of Schmitt et al. ( 2003b[ l. The 
reported discrepancy between their extrapolated ENLR sizes 
of luminous type 2 and type 1 QSOs is therefore obsolete. 

We find that the QSO continuum luminosity correlates much 
better with the ENLR size than does the total [Om] lumi- 
nosity, while the [Om] luminosity of the ENLR is totally 
uncorrelated with the QSO continuum luminosity. Thus, we 
argue that the [O m] luminosity is not a good proxy for the 
intrinsic luminosity of the QSO and that the amount of ion- 
ising photons does not determine the luminosity of an ENLR 
alone. Either the amount and distribution of gas is quite dif- 
ferent in the host galaxies or other factors such as the impact 
of radio jets play a additional role. This is supported by a 
trend that the ENLR luminosity increases with radio lumi- 
nosity. We therefore caution to use type 2 AGN to study the 
size-luminosity relation given that the [O m] luminosity is a 
secondary and less accurate proxy for the intrinsic luminos- 



ity of the nucleus (e.g. |Punsly & Zhang 2011 1. Addition 



ally, type 2 QSOs are often found to be radio-intermediate 
objects and have a much higher incidence of fiat-spectrum 
radio source (Lai & Ho 2010), which might cause a certain 
selection effect for type 2 QSOs if the radio jets enhance the 
[O m] luminosity relative to the bolometric AGN luminosity. 
Signatures of presumably rotational gas kinematics are 
identified in 6 host galaxies of which 4 are disc-dominated. 
Signatures of non-gravitational kinematics have been found 
in in 5 cases for specific regions within the EELRs among 
which 3 show radial velocity greater than >400kms _1 . 
Given the coincidences with extended radio sources, we 
argue that jet-cloud interactions are the most likely origin 
for the disturbances in the kinematics in the majority of 
cases. However, the detection rate and radial velocities 
of the outflows are much lower than reported for starburst 
galaxies or radio-loud AGN. 



One of the main results of our study is the low detection rate 
of gas outflows around low-redshift radio-quiet QSOs (RQQs), 
which is in agreement with similar results from type 2 QSOs and 
infrared faint Seyfert galaxies. Taken all together this is evidence 
that massive AGN-driven bulk outflows of gas are not ubiquitous 
around RQQs at low redshifts. Signatures of AGN feedback are 
mainly observed in radio-loud AGN which can expel the gas 
from the host galaxy or even drive large cavities in the intra- 
cluster medium. The lack of gas outflow signatures in RQQs 
may challenge some the theoretical AGN feedback scenarios in 
which the expulsion of gas is essential to stall star formation on 
short time scales. One explanation could be that the QSOs inves- 
tigated here, while among the most luminous at these low red- 
shifts, are still too faint compared to their more ultra-luminous 
high-redshift counterparts, to display the same strong signature 
of vigorous outflows in the ionised gas (Cano-Dfaz et al. 2012; 
Maio lino et al.|2012 i. But our RQQs may also be in late stages 
of their evolution, after most of the gas has either been heated 
or was already expelled from the host galaxies considering the 
presumed time delay between the peak of starburst and AGN ac- 
tivity of around a few 100 Myr ( |Li et al|2008j|Wild et al.|2010| i. 
On the other hand, the ENLR sizes, luminosities and kinemat- 
ics provide evidence for the general importance of radio jets to 



increase the covering factor of the ambient gas for the ionising 
photons of the nucleus in a pre-stage of feedback. To under- 
stand AGN feedback processes in RQQs therefore still remains 
an open issue. 
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Table 3. Broad-band imaging analysis results - Structural parameters and morphological host galaxy classification. 
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1.057 


1.068 


1 n/n /1 
1 .04 / ( 1 . /) 


1 //T /"2 AA 1 ,^ \ 

1.7 (3.96 kpc) 


0.40 


1 1'~ 




olJoo J 1033+2140 


CTACC 

oUbo \J ) 




1. lyj 


1 .002 




U . 2 (U.3 Kpc ) 


n Q/i 


C 1 

31 






SDSS (r') 




1.188 


1.140 




0'.'6(1.6kpc) 


0.82 


30° 




dp i inn i ciq 
ru 1 /UU+310 


rlSl Wrrtz (Tj4/M) 


(1) 


i onn 
t.8UU 


1. /4z 




1 n A /"A 1 

1.4 (0.2Kpc) 


n on 
0.8U 


1 1 

23 


M( V ) 


dp inn i aaa 
PCj 213U+UW 


UCT M7CDPT /CAAA^X7\ 

HS1 Wl'PCz (roUoW) 


(o) 


1.8/0 


1 QIC 

1.5/3 


1.833 (1.8) 


12 . (14 Kpc) 


n q a 
0.30 


An c 
00 


O 1 T"\ 




MTT C PiT7T ( U\ 




1 A 1 H 
1.41 / 


1 T70 
1. J /z 


1 3 An /"O /I A 

1.30U (2.4) 


1.5 (8.0 Kpc) 


n q 1 
0.31 


C 1 

3 1 




DUI 1 Ql 1 

r rlL loll 


UOT A PC /T7AO^A^7\ 

HS1 ACS (rozjW) 


(i) 


i.jyi 


1 CCO 

1. j3z 




1 . / (3.3 Kpc) 


n co 
0.32 


-70° 
11 


D(V) 




MTT C ( U\ 

IN 1 1 SUJ H 1 \ti ) 




l.ljo 


1 1 oo 

1. Izo 


1 nn 1 (c\ a\ 

i.oyi (y.o) 


. / (2.2 Kpc) 


n a 
0.48 


cno 

3y 






LK Apogee U4z \K) 




1.203 


1.248 


1.243 (8.3) 


1 . / (3.3 Kpc) 


n cn 
0.3y 


An c 
00 




HE2158-0107 f 


SDSS (z) 
NTT SOFI (H) 




0.884 


0.880 


0.880 (6.2) 


0'.'2 (0.7 kpc) 


0.93 


43° 


B? 


PG 2214+139 


NTT SOFI (H) 




1.014 


1.063 


0.950(1.8) 


2'.'9 (3.7 kpc) 


0.97 


64° 


B 




SDSS (;') 




1.825 


1.380 


1.299 (5.4) 


4'.'2 (5.4 kpc) 


0.94 


63° 






SDSS (r') 




2.326 


1.814 


1.684 (5.8) 


4'.'1 (5.3 kpc) 


0.94 


56° 




HE 2307-0254 


NTT SOFI (H) 




0.542 


0.542 


0.542 (6.2) 


0'.'2 (0.7 kpc) 


0.67 


35° 


B? 


PKS 2349-014 


HST WFPC2 (F606W) 


(8) 


2.649 


2.441 


2.436 (3.3) 


5'.'6 (16.5 kpc) 


0.81 


83' 


B/M 


HE 2353-0420 


LB Apogee U42 (tf ) 




0.940 


0.957 


0.939 (0.7) 


3'.'8 (13.9 kpc) 


0.49 


15° 


D(V) 



Notes. (a) References to the articles in which a morphological analysis of these data was initally presented. (b) Effective radius r c of the best-fitting 
realistic host-galaxy model given in apparent and physical scale. (c) Axis ratio of the minor to major axis of the best-fitting model. (d) Position angle 
of the major axis of the best-fitting model oriented such that 0° is North and 90° is East. (c) Our host galaxy morphological classification into: B 
(bulge-dominated), D (disc-dominated), B+D (bulge+disc component), M (Merger or strongly interacting), B?/D? (only tentative classifications), 
? (unconstrained). We append (V) to the classification, if it is due to a visual inspection of the image. (f) unresolved host galaxy 

Refer ences. (l)|Veilleux et al.| p009)>; (2) |Bahcall etaT1p996); (3)|McLure et al.| < fT999> ; (4) |Dunlop et~alj ( |2003) ; (5) |Bahcall etaL] (T997) >; 
(6) |Kim et alH2008> ; (7) |Jenkins et al.| ( |2005| l; (8) |Bahcal et al.|p995) 
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Table 4. PMAS IFU observational log. 



Object 


Obs. Date 


Exp. Time [s] 


Grating 


Wavelength Range (A) 


Seeing 


Obs. Condition 


PG 0026+129 


2002-09-05 


1x1200, 1x2400 


V300 


4600-7800 


l'.'l 


photometric 


IZw 1 


2002-09-03 


1x900, 2x1800 


V300 


3900-7200 


1'.'0-1'.'4 


mainly clear 


PG 0052+251 


2002-09-06 


2x1800 


V300 


4600-7800 


l'.'O 


photometric 


SDSS J0057+1446 


2008-09-05 


5x1800 


V600 


5100-6600 


175 


part, cloudy 


RBS 219 


2002-09-02 


1x1200, 1x1800 


V300 


4600-7800 


1'.'2-1'.'8 


photometric 


SDSS J0155-0857 


2008-09-07 


4x1800 


V600 


5100-6600 


1'.'4-1'.'8 


photometric 






2x1800 






172 


photometric 


HE 0157-0406 


2002-09-02 


1x1200, 1x1800 


V300 


5400-8600 


172-173 


photometric 


Mrk 1014 


2002-09-07 


1x1200, 1x1800 


V300 


5400-8600 


172-173 


photometric 




2008-09-06 


4x1800 


V600 


5100-6600 


174 


mainly clear 


SDSS J0836+4426 


2008-03-16 


3x1800 


V600 


5500-7000 


171 


part, cloudy 


SDSS J0948+4335 


2008-03-13 


3x1800 


V600 


5250-6750 


171 


cloudy 


SDSS Jl 13 1+2623 


2008-03-13 


3x1800 


V600 


5250-6750 


171 


part, cloudy 


PKS 1217+023 


2003-05-02 


4x1800 


V300 


5600-8900 


170-172 


mainly clear 


3C 273 


2003-05-03 


3x1200 


V300 


4350-7650 


174-179 


mainly clear 


HE 1228+0131 


2003-05-04 


4x1800 


V300 


4350-7650 


170-173 


part, cloudy 


SDSS J1230+6621 


2008-03-16 


2x1800 


V600 


5000-6500 


170 


mainly clear 


SDSS J1230+1 100 


2008-03-13 


1x1800 


V600 


5250-6750 


171 


cloudy 


PG 1427+490 


2003-05-01 


3x1200 


V300 


5600-8900 


079-172 


mainly clear 


SDSS J1444+0633 


2008-03-16 


1x1800 


V600 


5250-6750 


173-175 


cloudy 


HE 1453-0303 


2003-05-02 


2x1800 


V300 


5600-8900 


172-175 


mainly clear 




2003-05-03 


2x1800 






173-176 


mainyl clear 


PKS 1545+21 


2003-05-03 


3x1800 


V300 


5600-8900 


172-175 


mainly clear 


PG 1612+261 


2008-09-07 


3x1800 


V600 


5100-6600 


171 


mainly clear 


SDSS J1655+2146 


2008-03-16 


3x1200 


V600 


5000-6500 


173-176 


part, cloudy 






1x1800 






170 


photometric 


PG 1700+518 


2003-05-02 


2x1800 


V300 


5600-8900 


079-172 


photometric 




2003-05-03 


2x1800 






176 


mainly clear 


PG 2130+099 


2002-09-04 


2x2400 


V300 


3900-7200 


173 


photometric 


PHL 1811 


2002-09-06 


2x1800 


V300 


5400-8600 


172 


photometric 


HE 2158-0107 


2002-09-06 


2x1800 


V300 


5400-8600 


171 


photometric 


RBS 1812 


2002-09-04 


1x2400 


V300 


4600-7800 


175 


photometric 




2002-09-05 


2x2400 






176-270 


mainly clear 


PG 2214+139 


2002-09-02 


2x2400, 1x1800 


V300 


3900-7200 


173-174 


photometric 


HE 2307-0254 


2002-09-06 


2x1800 


V300 


5400-8600 


172 


photometric 


PKS 2349-014 


2002-09-06 


3x900 


V300 


5400-8600 


170-171 


photometric 


HE 2353-0420 


2002-09-07 


2x1800 


V300 


5400-8600 


171-174 


photometric 
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Table 7. Kinematics and line fluxes for specific regions of the EELR. 



Object 11 


Av r b 


C line° 


w 


[O m] ,15007 Ha 


[N n] /165 85 Classification 




[kms- 1 ] 


[kms- 1 ] 






10 16 ergs 1 cm 2 





PG 0026+129 


A 


48 ± 34 


158 ± 14 


2.4 ± 0.9 


25.1 ± 1.0 






AGN 




B 


136 + 32 


153 + 9 


2.5 + 0.4 


19.8 + 0.5 






AGN 


PG 0050+124 


A 


-312 + 27 


83 + 24 






9.3+ 1.4 


7.9 + 0.3 


Hn 




B 


-475 ± 28 


18 + 23 






6.9 + 2.8 


7.7 ± 0.4 


Hn 




C 


-215 + 36 


80 + 30 






9.8+1.9 


5.8 + 0.7 


Hn 


PG 0052+251 


A 


-127 + 36 


110 + 31 


3.0 + 0.5 


9.4 ± 0.2 


9.5 + 0.4 


< 2.5 


Hn 




B 


-232 ± 32 


82 + 27 


6.3+1.1 


38.5 ± 0.2 


22.1 + 0.8 


< 3.4 


Hn? 




C 


-264 + 41 


65 + 28 


5.5 + 1.2 


7.0 + 0.4 


14.1 ± 1.1 


< 5.3 


AGN 


RBS 219 


A 


-95 ± 40 


76 + 22 






12.1 ± 1.2 


< 2.5 


Hn 




B 


17 + 39 


0+25 






5.8+1.7 


< 4.6 


Hn 




C 


-105 + 47 


0+39 






4.7 ± 2.9 


< 8.3 


Hn 


SDSS J0155-0857 


A 


134 ± 40 


73 ± 12 


0.4 ± 0.3 


3.5 ± 0.3 






AGN 




B 


-7 + 42 


129+ 13 


1.7 + 0.8 


7.0 + 0.6 






AGN 


Mrk 1014 


A 


- 1 00 + 28 


155 ± 22 


9.1 ± 1.8 


53.9 ± 0.5 


26.1 ± 5.1 d 


14.6 ± 0.7 


AGN 




B 


67 ± 32 


229 ± 17 


14.0 + 3.5 


40.1 ±0.9 


40.2 ± 10.0 d 


14.2 ± 0.8 


AGN 




C 


-222 + 31 


113 + 28 


5.0+1.0 


23.0 ± 0.5 


14.2±2.7 d 


6.8+1.0 


AGN 


SDSS J0846+4426 


A 


-88 ± 46 


132 + 3 


18.2 + 2.2 


198.2 + 2.3 






AGN 




B 


-229 ± 46 


143 + 3 


14.2+1.4 


146.0+ 1.5 






AGN 




C 


-154 + 50 


130+8 


< 5.4 


58.6+1.8 






AGN 


SDSS J1131+2623 


A 


-44 ± 69 


117 + 3 


1.9 + 0.9 


32.6 ± 0.6 






AGN 




B 


-32 ± 68 


87 + 4 


1.6 + 0.3 


25.7 ± 0.4 






AGN 




C 


-424 ± 72 


54+8 


< 1.3 


8.6 + 0.4 






AGN 




D 


131 + 76 


105 ± 10 


< 2.4 


13.9 + 0.8 






AGN 


PKS 1217+023 


A 


45 + 35 


75 + 27 


2.5 ± 0.6 


9.3 ± 0.3 






AGN 




B 


-61 ± 37 


80 + 27 


3.4+1.1 


19.4 + 0.7 


6.1+0.8 


2.8 ± 0.5 


AGN 


SDSS Tl ?3f)+(Sfi? 1 


A 


299 + 24 


U4 + 5 


4 n + n o 


3 + Q 
jj.j in \j.y 






AGN 




B 


67 ± 23 


80 + 5 


2.9 + 0.4 


19.5 + 0.5 






AGN 


PKS 1545+210 


A 


309 ± 47 


93 + 29 


4.5 ± 2.4 


36.9 ± 1.8 






AGN 




B 


-134 + 64 


212 + 51 


6.0 + 4.1 


35.5+4.7 






AGN 


PG 1612+261 


A 


70 + 29 


77 + 5 


< 18.9 


111.0 + 4.7 






AGN 




B 


54 + 28 


91 + 3 


< 8.2 


289.1 +2.0 






AGN 




C 


248 ± 28 


82 + 4 


18.2 + 3.7 


199.8+1.2 






AGN 




D 


-92 + 33 


105 + 8 


< 4.8 


29.3 ± 1.2 






AGN 


SDSS J1655+2146 


A 


-85 ± 26 


161 + 2 


1.3 + 0.8 


127.9 + 0.9 






AGN 




B 


-85 + 27 


183 + 3 


20.4 ± 2.6 


319.3 + 3.0 






AGN 




C 


-78 + 31 


197 + 7 


8.1 ± 1.0 


44.3 ± 1.2 






AGN 




D 


-128 + 30 


94 + 5 


< 4.1 


43.5 ± 1.4 






AGN 


PG 1700+518 


A 


-889 ± 67 


158 + 39 


2.1+0.6 


9.2 + 0.8 






AGN 


PG 2130+099 


A 


-99 ± 68 


180+105 


5.9 + 3.0 


3.5 ± 0.7 


7.7+1.6 


3.8 ± 0.4 


Hn 




B 


-262 ± 36 


176 + 32 


10.1 + 2.8 


22.1 ± 1.3 


21.6+1.1 


9.5 ± 0.3 


AGN? 


PHL 1811 


A 


49 + 27 


33+ 15 


3.4 + 0.4 


3.4 ± 0.4 


14.6+ 1.0 


5.7 ± 0.5 


Hn 


HE 2158-0107 


A 


42 + 44 


69 + 26 


8.1 ± 1.4 


39.6+ 1.5 


11.5+ 1.6 


< 3.8 


AGN 




B 


-50 + 39 


39 + 23 


11.2+1.8 


108.5 ± 1.9 


28.4+3.6 


< 8.2 


AGN 




C 


-50 ± 47 


0+30 


< 3.6 


18.7+1.2 






AGN 


PKS 2349-014 


A 


-33 ± 49 


172 ± 24 


4.2 + 0.8 


6.9 + 0.3 


10.3 ± 0.7 


5.7 ± 0.4 


Hn 




B 


-269 ± 57 


299 + 21 


5.6 + 2.1 


20.8 + 1.1 


11.7 + 2.2 


13.2+ 1.2 


AGN 




C 


97 + 39 


146 ± 22 


12.2+1.2 


96.6 ± 0.9 


34.0 ± 2.2 


36.8+ 1.6 


AGN 




D 


418 + 47 


+ 26 


1.3 + 0.1 


5.8 + 0.3 






AGN 


HE 2353-0420 


A 


85 + 51 


225 ± 17 


7.6+1.5 


77.9 ± 2.0 


22.5+ 1.6 


18.0+ 1.0 


AGN 




B 


66 + 46 


189+18 


6.0+1.3 


116.5 + 1.8 


32.5 ± 2.2 


26.6+ 1.6 


AGN 



Notes. <a) Name of QSOs followed with the alphabetically ordered character of the specific co-added spectra. <b) Radial velocity offset with respect 
to the rest-frame of the nuclear [O m] line. Negative values corresponds to a blue shift. <c) Emission line dispersion after the instrumental dispersion 
was subtracted in quadrature. A zero dispersion means that the line was unresolved. (d) The Ha line was strongly affected by telluric absorption 
and its flux was estimates from the H/3 line assuming zero dust extinction (Ha/H/3=2.87). 



